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We report on the activities of the "QCD Tools for heavy flavors and new physics searches" working group of the Run II 
Workshop on QCD and Weak Bosons. The contributions cover the topics of improved parton showering and comparisons 
of Monte Carlo programs and resummation calculations, recent developments in PYTHIA, the methodology of measuring 
backgrounds to new physics searches, variable flavor number schemes for heavy quark electro-production, the underlying 
event in hard scattering processes, and the Monte Carlo MCFM for NLO processes. 



1. Overview 

The task of the "QCD Tools for heavy flavors and 
new physics searches working group" was to evaluate 
the status of the tools - invariably computer programs 
that simulate physics processes at colliders - that 
are being used to estimate signal and background 
rates at the Tevatron, and to isolate areas of con- 
cern. The contributions presented here cover several 
topics related to that endeavor. It is hoped that 
the next period of data-taking at the Tevatron will 
reveal indirect or direct evidence of physics beyond 
the Standard Model. The precise measurement of the 
W boson mass Mw and its correlation with the top 
quark mass rrit is one example of an indirect probe 
of the Standard Model. The production of a light 
Higgs boson in association with a. W or Z boson 
is an obvious example of a direct one. While both 
measurements are related to electroweak symmetry 
breaking, it requires a quantitative understanding of 
perturbative and non-perturbative QCD to interpret 
data. 

Because of the importance of the Mw measurement, 
and since gauge boson production in association with 
jets is a serious background in many new physics 
searches, much effort was devoted to understanding 
gauge boson production processes. It is well known 
that the emission of many soft gluons has a profound 
effect on the kinematics of gauge boson production. 
Two calculational methods have been used to com- 
pare "theory" with data: (1) analytic resummation 
of several series of important logarithms, and (2) 
parton showering based on DGLAP-evolved parton 
distribution functions. Here, there are reports on our 
understanding of both, and improvements. Note also 
that diboson production is often a background too. 

In the Standard Model, and its minimal super- 
symmetric extension, the mechanism that generates 



mass for the electroweak gauge bosons also generates 
fermion mass. From an agnostic point of view, the 
fact that the W and Z bosons and the top quark 
have roughly similar masses, and these masses are 
quite disparate from, say, the electron or neutrino 
masses, is some evidence that heavy flavor is related 
to electroweak symmetry breaking. Many of the search 
strategies for Run II rely on tagging c and b quarks or t 
leptons. For this reason, there are several contributions 
regarding issues of determining backgrounds in Run II. 

2. Performing parton showering at Next-to- 
Leading-Order Accuracy 

by S. Mrenna 

2.1. Introduction 

In the near future, experiments at the Tevatron 
will search for evidence of physics that supersedes the 
standard model. Important among the tools that will 
be used in these searches are showering event gener- 
ators or showering Monte Carlos (SMC's). Among 
the most versatile and popular of these are the Monte 
Carlos HERWIG[§, ISAJET[|l, and PYTHIA[|. 
SMC's are useful because they accurately describe the 
emission of multiple soft gluons, which is, in effect, 
an all orders problem in QCD. However, they only 
predict total cross sections to a leading order accuracy, 
and, thus, can demonstrate a sizeable dependence on 
the choice of scale used for the parton distribution 
functions (PDF's) or coupling constants (particularly 
as). Also, in general, they do not translate smoothly 
into kinematic configurations where only one, hard 
parton is emitted. In distinction to SMC's are certain 
analytic calculations which account for multiple soft 
gluon emission and higher order corrections to the hard 
scattering. These resummation calculations, however, 
integrate out the kinematics of the soft gluons, and, 
thus, are limited in their predictive power. They 
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can, for example, describe the kinematics of a heavy 
gauge boson produced in hadron coUision, but can- 
not predict the number or distribution of jets that 
accompany it. However, searches for new physics, 
either directly or indirectly through measurements 
of precision electroweak observables, often demand 
detailed knowledge of kinematic distributions and jet 
activity. Furthermore, W^+jets (and Z+jets) processes 
are often backgrounds to SUSY or technicolor signa- 
tures, and we demand a reliable prediction of their 
properties. Here, we report on recent progress in 
improving the predictive power of showering Monte 
Carlos by incorporating the positive features of the 
analytic resummation calculations into the showering 
algorithms. In the ensuing discussion, we focus on 
the specific example of W boson production at a 
hadron collider, when the W decays leptonically. The 
results apply equally well toj*,Z and Higgs bosons (or 
any heavy, color-singlet particle) produced in hadron 
collisions. 

2.2. Parton Showers 

SMC's are based on the factorization theorem [|^, 
which, roughly, states that physical observables in 
QCD are the product of short-distance functions and 
long-distance functions. The short-distance functions 
are calculable in perturbation theory. The long- 
distance functions are fit at a scale, but their evolution 
to any other scale is also calculable in perturbation 
theory. 

A standard application of the factorization theorem 
is to describe W boson production at a pp collider 
at a fixed order in a^. The production cross section 
is obtained by convoluting the partonic subprocesses 
evaluated at the scale Q with the PDF's evaluated at 
Q. The partons involved in the hard collision must 
be sufficiently virtual to be resolved inside the proton, 
and a natural choice for the scale Q is Q — Mw [ 
|]. However, the valence quarks in the proton have 
virtualities at a much lower scale Qq of the order 
of 1 GeV. The connection between the partons at 
the low scale Qo and those at the high scale Q is 
described by the DGLAP evolution equations [ |^. 
The DGLAP equations include the most important 
kinematic configurations of the splittings a —^ be, 
where a, b and c represent different types of partons in 
the hadron {q,g, etc.). Starting from a measurement 
of the PDF's at a low scale Qq, a solution of the 
DGLAP equations yields the PDF's at the hard scale 
Q. Equivalently, starting with a parton c involved 
in a hard collision, it is also possible to determine 
probabilistically which splittings generated c. In the 
process of evolving parton c back to the valence quarks 
in the proton, a number of spectator partons (e.g. 



parton b in the branching a — > be) are resolved. These 
partons constitute a shower of soft and/or coUinear 
jets that accompany the VF-boson, and influence its 
kinematics. 

The shower described above occurs with unit prob- 
ability and does not change the total cross section for 
W^-boson production calculated at the scale Q [ 0]. 
The showering can be attached to the hard-scattering 
process based on a probability distribution after the 
hard scattering has been selected. Once kinematic cuts 
are applied, the transverse momentum and rapidity 
of the Tl^-boson populate regions never accessed by 
the differential partonic cross section calculated at a 
fixed order. This is consistent, since the fixed-order 
calculation was inclusive (i.e., pp ^ W + X) and was 
never intended to describe the detailed kinematics of 
the VF-boson in isolation. The parton shower, in effect, 
resolves the structure of the inclusive state of partons 
denoted as X. In practice, the fixed order partonic 
cross section (without showering) can still be used to 
describe properties of the decay leptons as long as the 
observable is well defined (e.g., the number of leptons 
with central rapidity and high transverse momentum, 
but not the distribution of transverse momentum of 
the W). 

Here, we focus on the case of initial state gluon 
radiation. More details can be found in Ref. for 
example. Showering of the parton b with momentum 
fraction x resolved at the scale = e* is driven by a 
Sudakov form factor exp(— S"), such as[ ^ 

x' = x/z, (1) 

which is implemented in PYTHIA, and the formally 
equivalent expression [ [lO|] 

Mt') fa{x,t) 

h{x,t') A{t) ' 

A(t') = exp|^-^* J'~' dt"dz^^^^Pa^,,{z)^ 

which is implemented in HERWIG. In the above 
expressions, to is a cutoff scale for the showering, P 
is a DGLAP splitting function, and fi is a parton 
distribution function. The Sudakov form factor pre- 
sented here is a solution of the DGLAP equation, 
and gives the probability of evolving from the scale 
= e* to Q'^ = e* with no resolvable branching. 
The Sudakov form factor contains all the information 
necessary to reconstruct a shower, since it encodes 
the change in virtuality of a parton until a resolvable 
showering occurs. A parton shower is then an iterative 
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solution of the equation r — exp(— S'), where r is a 
random number uniformly distributed in the interval 
[0, 1], until a solution for Q' is found which is below a 
cutoff. For consistency, the cutoff should represent the 
lowest scale of resolvable emission Qo- The evolution 
proceeds backwards from a large, negative scale — IQ'^I 
to a small, negative cutoff scale — |Qol- 

After choosing the change in virtuality, a particu- 
lar backwards branching and the splitting variable z 
are selected from the probability function based on 
their relative weights (a summation over all possible 
branchings a — » 6c is implied in these expressions). 
The details of how a full shower is reconstructed in the 
PYTHIA Monte Carlo, for example, can be found in 
Ref. [ 1^ . The structure of the shower can be complex: 
the transverse momentum of the VF-boson is built up 
from a whole series of splittings and boosts, and is 
known only at the end of the shower, after the final 
boost. 

The SMC formulation outlined above is fairly in- 
dependent of the hard scattering process considered. 
Only the initial choice of partons and possibly the high 
scale differs. Therefore, this formalism can be applied 
universally to many different scattering problems. In 
effect, soft and collinear gluons are not sensitive to the 
specifics of the hard scattering, only the color charge 
of the incoming partons. 

2.3. Analytic Resummation 

At hadron colliders, the partonic cross sections 
can receive substantial corrections at higher orders 
in as- This affects not only the total production 
rate, but also the kinematics of the W boson. At 
leading order (a^), the M^-boson has a S{Q'^) dis- 
tribution in Q^. At next-to-leading order, the 
real emission of a single gluon generates a contri- 
bution to da/dQ^ that behaves as Q^^as{Qx) and 
ln((3^/Qy) while the leading order, soft, 
and virtual corrections are proportional to —d{Q^). 
At higher orders, the most singular terms follow the 

pattern of E™Jo Ih^IQVQt) = = 

y". The logarithms arise from the incomplete can- 
cellation of the virtual and real QCD corrections, but 
this cancellation becomes complete for the integrated 
spectrum, where the real gluon can become arbitrarily 
soft and/or collinear to other partons. The pattern 
of singular terms suggest that perturbation theory 
should be performed in powers of V" instead of a". 
This reorganization of the perturbative series is called 
resummation. 

The first studies of soft gluon emission resummed the 
leading logarithms [ 11, leading to a suppression 
of the cross section at small Qt- The suppression 
underlies the importance of including sub-leading log- 



arithms [ The most rigorous approach to the 

problem of multiple gluon emission is the CoUins- 
Soper-Sterman (CSS) formalism for transverse mo- 
mentum resummation [ which resums all of the 
important logarithms. This is achieved after a Fourier 
transformation with respect to Qt in the variable b, so 
that the series involving the delta function and terms 
V" simplifies to the form of an exponential. Hence, 
the soft gluon emission is resummed or exponentiated 
in this 6-space formalism. Despite the successes of 
the 6-space formalism, there are drawbacks: the soft 
gluon dynamics are integrated out, and the Sudakov 
form factor is a Fourier transform. 

The CSS formalism was used by its authors to 
predict both the total cross section to NLO and the 
kinematic distributions of the VF-boson to all orders 
[ |l^ at hadron colliders. A similar treatment was 
presented using the AEGM formalism [ ^ , that does 
not involve a Fourier transform, but is evaluated 
directly in transverse momentum Qt space. When 
evaluated at NLO, the two formalisms are equivalent 
to NNNL order in as, and agree with the fixed order 
calculation of the total cross section [ A more 

detailed numerical comparison of the two predictions 
can be found in Ref. [ ^ . 

Recently, the AEGM formalism has been re- 
investigated, and an approximation to the 6-space 
formalism has been developed in Qr-space which 
retains its predictive features [ |l^ (see also the recent 
eprint [ ^). This formulation does have a simple, 
physical interpretation, and can be used to develop 
an alternate algorithm for parton showering which in- 
cludes higher-order corrections to the hard scattering. 
For this reason, we focus on the QT~space formalism. 
To NNNL accuracy, the Qt space expression agrees 
exactly with the 6-space expression, and has the form 



da{hih2 y(*)x) 

dg2 dQ^ dy 



-W{Qt-,Q, Xi,X2) 



dQ^ 

+Y{Qt,Q,xi,X2). (3) 



In this expression, Q, Qt and y describe the kinematics 
of the boson V, the function Y is regular as Qt 
and correjrts for the soft gluon approximation, and the 
function W has the form: 



(c®/)(xi,QT)(c®/)(a;2,( 



(4) 



where 



SiQT,Q) = 
dp^ 

Ql 



(5) 
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and 



{Cji® fi/h,) {xi,ti) = 



/ fic,Kfi,QT)M.(a,QT) 



(6) 



_ff is a function that describes the hard scattering, and 
A, B, and C are calculated perturbatively in powers 
of a,: 



n=0 



{A,B,C) 



(n) 



(the first non-zero terms in the expansion of A and 
B are for n = 1). The functions C'-"-' are the Wilson 
coefficients, and are responsible for the change in the 
total production cross section at higher orders. In fact, 
(C ® /) is simply a redefinition of the parton distri- 
bution function obtained by convoluting the standard 
ones with an ultraviolet-safe function. 

Ignoring Y and other kinematical dependence, 
Eq. can be rewritten as: 

da{hih2 -> WX) ( d 



dQl 



^-SiQT,Q) ^ 



where 



R 



and 



{C(g>f) ixi,QT)iC(g>f) (x2,Qr) 
iC(S>f) ixi,Q) {C®f) {X2,Q) 



dxi 



(7) 



(8) 



ai=K I ^ (C ® /) ixi,Q) (C ® /) ix2, Q). 

Xl 

The factor cti is the total cross section to a fixed order, 
while the rest of the function yields the probability that 
the W^-boson has a transverse momentum Qt- 

At leading order, the expression for the production 
of an on-shell VF-boson simplifies considerably to: 



da{hih2 -> WX) 



dQl 



dQl 



f{x^,Q)f{x2,Q) 



(9) 



2.4. A modified showering algorithm 

The primary result of this analysis is to exploit 
the expression for the differential cross section, which 
has the form of a leading order cross section times a 
backwards evolution, to incorporate NLO corrections 
to the parton shower. We generalize the function 
A(t) / f{x,t) X f{x, t') /A{t') of the standard backwards 

showering algorithm to \/w (the square root appears 
because we are considering the evolution of each parton 
line individually). 

To implement this modification in a numerical pro- 
gram, like PYTHIA, we need to provide the new, mod- 
ified PDF (mPDF) based on the Wilson coefficients. 
At leading order, the only Wilson coefficient is C^j^ = 
SijS(l — z), and we reproduce exactly the standard 
showering formulation. For W^-boson production at 
NLO, the Wilson coefficients C are: 



C!jl} = 6,,\1{1-z) + 1{7t'-8)6{1-z)}, (10) 



a 



(1) 



1 



J9 2' 

To NLO, the convolution integrals become: 



as{^l) 1 2 



(11) 



{C^ fi){x,fi) ^ f,ix,fi) 1 + 



TT 3 



c^si^J■) f dz 
z 



Kl - z)fi{x/z,n) 



+ 2(1 - z)fg{x/z,^i) 



and fg{x,ii) is unchanged. The first term gives the 
contribution of an unevolved parton to the hard scat- 
tering, while the other two contain contributions from 
quarks and gluons with higher momentum fractions 
that split q qg and g qq, respectively. 

We are assuming that the Sudakov form factor 
used in the analytic expressions and in the SMC are 
equivalent. In fact, the integration over the quark 
splitting function in A{Q) yields an expression similar 
to the analytic Sudakov: 



CTO 



dxi 

Xl 



f{xi,Q)f{x2,Q), 



where k contains physical constants and we ignore 
the function Y for now. The expression contains 
two factors, the total cross section at leading order 
ctq, and a cumulative probability function in that 
describes the transverse momentum of the VF-boson. 
The term e~^^^f{x,QT)/fix,Q) is of the same form 
as the Sudakov form factor of Eq. (2) and, hence, to 
that of Eq. (1). 



dzCp 



Cf In 



l + z' 
1- z 

2 



3/2(1 - 2Zrn) 



where 



AWln(QVQT)+S^'\ (13) 



is an infrared cutoff, terms 



(Q+Qt) 

of order z„i and higher are neglected, and the z 
dependence of the running coupling has been ignored 
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[ H. Note that the coefficients A^^) (Cf) and B^^) 
(— 3/2Cf) are universal to qq annihilation into a color 
singlet object, just as the showering Sudakov form 
factor only knows about the partons and not the 
details of the hard scattering. For gg fusion, only the 
coefficient A^^^ (3) is universal. In general, at higher 
orders, the analytic Sudakov is sensitive to the exact 
hard scattering process. 

While the Sudakov form factors are similar, there is 
no one-to-one correspondence. First, the Qr-space 
Sudakov form factor is expressed directly in terms 
of the Qt of the heavy boson, while, in the SMC's, 
the final Qt is built up from a series of branchings. 
Secondly, the integral on the left of Eq. (13) is positive 
(provided that Zm < \), while the analytic expression 
on the right can become negative. This is disturbing, 
since it means sub-leading logarithms (proportional 
to B) are dominating leading ones. In the exact 
SMC Sudakov, the kinematic constraints guarantee 
that A((5) < 1. In this sense, the Sudakov in the 
SMC is a more exact implementation of the analytic 
one. Nonetheless, the agreement apparent between the 
analytic and parton shower expressions is compelling 
enough to proceed assuming the two Sudakov form 
factors are equivalent. 



2.5. Hard Emission Corrections 

The SMC and resummation formalisms are opti- 
mized to deal with kinematic configurations that have 
logarithmic enhancements L. For large Qt — Q, 
there are no such enhancements, and a fixed order 
calculation yields the most accurate predictions. The 
region of medium Qt, however, is not suited to either 
particular expansion, in a'^L'^ or a". 

The problem becomes acute for SMC's. In the 
standard implementation of SMC's, the highest Qt 
is set by the maximum virtuality allowed, Q — Myy 
in our example, so that the region Qt > Q is never 
accessed. However, at Qt > Q, the fixed order calcula- 
tion is preferred and yields a non-zero result, so there 
is a discontinuity between the two predictions. This 
behavior does not occur in the analytic resummation 
calculations, because contributions to the cross section 
that are not logarithmically enhanced as Qt — > 
are added back order-by-order in a^. This procedure 
corrects for the approximations made in deriving the 
exponentiation of soft giuon emission. The correction 
is denoted Y. If the coefficients A and B are calculated 
to high-enough accuracy, one sees a relatively smooth 
transition between Eq. ^ and the NLO prediction at 
Qt = Q- In the Qr^spa-ce calculation, this matching 
between the two calculations at Qt = Q is guaranteed 



at any order. The function Y has the form 



Y{Qt,Q,Xi,X2) = — ~^2^ 

faiiu Q) Ril\QT, Q,^,^) /b(6, Q). (14) 

For W OT Z boson production, the a — q,b — q 
component of R at first order in as is 



Rq\^ — Cp — 



Q'' 



r + [u - Q') 



2^2 



iu 



■S{s + i+u-Q^) 



-Pg^g{zB)Sil- ZA)-iA^ B). (15) 



The invariants s,t and u are defined in terms of 
z,Q,Qt- 



t/Q^ = l- l/zBx 1 + Q^/Q^ 



ii/Q^ = 1 - l/z.4^1 + Q|/g2. 



(16) 



The term in R proportional to the delta function 
is simply the squared matrix element for the hard 
emission, while the terms proportional to Q^^ are the 
asymptotic pieces from W. 

We would like to include similar corrections into the 
SMC. However, this is not entirely straightforward. 
Though it is not obvious from Eq. (14), the {a = g,b ^ 
g-t-permutations) components are negative for Qt < 
Q, though the sum Y is positive. Retaining negative 
weights in an intermediate part of the calculation 
is not a problem in principle. We can artificially 
force the negative weights to be positive, and then 
include the correct sign of the weight when filling 
histograms, for example. However, this would involve 
some modification to the PYTHIA code used in this 
study. 

A pragmatic approach is to ignore the negative 
weights entirely, and multiply the exact VF-|-parton 
cross sections by a factor so that their sum reproduces 
the Qt distribution and normalization of the analytic 
Y piece. For the Tevatron in Run I, we find that the 
multiplicative factor /cor — |(Qt/50)^x (H-(5t/25) 
reproduces the correct behavior for Qt < 50 GeV. 
For Qt > 50 GeV, the uncorrected W+ parton cross 
sections are employed. Since the matching between 
the "resummed" and "fixed order" calculations is now 
occurring at Qt — 50 GeV instead of Qt = Mw, 
we further limit the maximum virtuality of showering 
to 50 GeV. This is in accord with the fact that the 
"resummed" part of the analytic calculation becomes 
negative around Qt — 50 GeV. This choice does have 
some effect on the overall normalization of the parton 
showering component. 
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At this point, it is useful to compare the scheme 
outhned above to other approaches at improving the 
showering algorithm. One scheme is based on phase- 
space splitting of a NLO matrix element into a piece 
with LO kinematics and another with exclusive NLO 
kinematics [ ^ The separation depends on an 

adjustable parameter that splits the phase space. In 
the approach of Ref. [ |2^, the separation parameter 
is tuned so that the contribution with LO kinematics 
vanishes. The resultant showering of the term with ex- 
clusive NLO kinematics can generate emissions which 
are harder than the first "hard" emission, which is 
not consistent. More seriously, physical observables 
are sensitive to the exact choice of the separation 
parameter (see the discussion in Ref. [ ^ regarding 
Q^^). Furthermore, the separation parameter must be 
retuned for different processes and different colliders. 
This scheme is guaranteed to give the NLO cross 
section before cuts, but does not necessarily generate 
the correct kinematics. 

The other scheme is to modify the showering to 
reproduce the hard emission limit [ |2^. While 
this can be accomplished, it does so at the expense of 
transferring events from low Qt to high Qt- There 
is no attempt to predict the absolute event rate, 
but only to generate the correct event shapes. In 
some implementations, the high scale of the showering 
is increased to the maximum virtuality allowed by 
the collider energy. This is contrary to the analytic 
calculations, where the scale Q — M^y, for example, 
appears naturally (in the choice of constants Ci , C2 
and C3 which eliminate potentially large logarithms). 
This scheme will generate the correct hard limit, but 
will not generate the correct cross section in the soft 
limit. 

2.6. Numerical results 

For our numerical results, we predict the Qt distri- 
bution of W and Z bosons produced at the Tevatron 
in Run I. The modified PDF (mPDF) was calculated 
using CTEQ4M PDF's. These distributions are in 
good agreement with analytic calculations, but the 
shape and overall normalization cannot be predicted 
accurately by the standard showering algorithm. Some 
of the alternative showering algorithms reproduce the 
shape, but not the overall normalization. Secondly, we 
discuss jet properties for the same processes, which 
are not significantly altered from the predictions of 
the standard showering algorithm. These cannot be 
predicted by analytic calculations. 

In Fig. 0(a), the transverse momentum of the 
W boson (solid line) as predicted by the algorithm 
outlined above is shown in comparison to D0 data[ 
p7| (crosses). The theoretical distribution has been 



passed through the CMS detector simulation. [| As 
in analytic calculations, the position of the peak in 
the Qt distribution from parton showering depends 
on non-perturbative physics [ In PYTHIA, 

this is implemented through a Gaussian smearing of 
the transverse momentum of the incoming partons. 
To generate this plot, we have changed the default 
Gaussian width from .44 GeV to 2.0 GeV, which 
is more in accord with other analyses. This is the 
value used in all subsequent results. Because of 
the necessity of reconstructing the missing Et in W 
boson decays to leptons, the smearing of the Qt 
distribution is significant, and the agreement between 
the prediction and data is not a rigorous test of the 
modified showering algorithm. Fig. ^b)shows the 
comparison of the CDF Drell-Yan data [ Q near the 
Z° peak to the modified showering prediction. While 
there is a problem with the overall normalization, the 
shape agreement is very good. We note that there is 
also a problem with the overall normalization of the 
analytic resummation predictions]^ 

Given all the effort necessary to improve the shower- 
ing, it is reasonable to ask if the similar results would 
have been obtained by simply renormalizing the usual 
predictions to the NLO rate, i.e. using PYTHIA but 
applying a constant if-factor at the end. In W boson 
production, the relative size of the Qt distributions 
vary by as much as 10% in the important regions 
of small and medium Qt- Of course, the effect is 
much larger for the large Qt region where there is 
almost no rate from the standard parton showering. 
If one is worried about precision measurements or 
is applying kinematic cuts that bias the large Qt 
region, then standard parton showering can yield 
misleading results. In most cases, however, it appears 
to be perfectly reasonable to renormalize the parton 
showering results to the total NLO cross section. We 
have also checked if our new showering algorithm has 
an impact on jet properties. For W and Z boson 
production, there are only minor differences, which is 
expected since the Wilson coefficients for W and Z 
boson production are nearly unity. In general, we do 
not expect any major changes from using the modified 
PDF's, since the showering depends on the ratio of 
the modified PDF's evaluated at two different scales, 
which is not as sensitive to the overall normalization 
of the PDF. 



^Special thanks go to Cecilia Gerber, for making the code 
portable, and to Michael Seymour for explaining how to properly 
use it. 

■^Csaba Balazs, private communication. 
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W Boson Production at the Tevatron, Run I 



Z Boson Production at tlie Tevatron, Run I 




Figure 1. (a) The prediction of the W boson transverse momentum distribution in Run I at the Tevatron (solid 
hne) compared to the D0 data. The prediction includes the effects of the modified parton distribution functions, 
the correction to the hard scattering process, and a primordial kx of 2.0 GeV; (b) The prediction of the Z boson 
transverse momentum distribution in Run I at the Tevatron (solid line) compared to the CDF data. The prediction 
includes the effects of the modified parton distribution functions, the correction to the hard scattering process, and 
a primordial kx of 2.0 GeV. 



2.7. Conclusions 

We have presented a modified, parton showering 
algorithm that produces the total cross section and 
the event shapes beyond the leading order. These 
modifications are based on the Qj^-space resummation. 
The parton showering itself is modified by using a new 
PDF (called mPDF) which encodes some information 
about the hard scattering process. Simultaneously, 
the explicit, hard emission is included, but only after 
subtracting out the contribution already generated 
by the showering: this correction is called Y. The 
presence of Y yields a smooth transition from the 
parton showering to single, hard emission. We mod- 
ified the PYTHIA Monte Carlo to account for these 
corrections, and presented comparisons with Run I W 
and Z boson data. 

The scheme works well for the cases considered in 
this study, and the correct cross sections, transverse 
momentum distributions, and jet properties are gener- 
ated. We have compared our kinematic distributions 
to the case when the results of the standard showering 
are multiplied by a constant iC-factor to reproduce the 
NLO cross section. We find variations on the order of 
10% for small and medium transverse momentum. 



There are several effects which still need study. We 
have not included the exact distributions for the decay 
of the leptons [ ^ for W and Z production, which 
are resummcd differently. It is straightforward to 
include such effects. In the theoretical discussion and 
numerical results, we have focussed on initial state 
radiation, but our results should apply equally well for 
final state radiation. The situation is certainly simpler, 
since final state radiation does not require detailed 
knowledge of the fragmentation functions. Also, the 
case when color flows from the initial state to the final 
state requires study. A resummed calculation already 
exists for the case of deep inelastic scattering [ 
and much theoretical progress has been made for heavy 
quark production 
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We believe that the modified 
showering scheme outlined in this study generalizes 
beyond NLO, just as the analytic calculations can be 
calculated to any given order. For example, we could 
include hard 14^ -I- 2 jet corrections [ ^ to Y. For 
consistency, however, higher order terms {A and B) 
may also need to be included in the Sudakov form 
factor. 

The modified PYTHIA subroutines used in this 
study and an explanation of how to use them are 
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available at the following URL: 

moose . ucdavis . edu/mrenna/ shower . html. 
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3. Recent Progress in PYTHIA 
by T. Sjostrand 

3.1. Introduction 

A general-purpose generator in high-energy physics 
should address a number of physics aspects, such as: 

• the matrix elements for a multitude of hard 
subprocesses of interest, 

• the convolution with parton distributions to 
obtain the hard-scattering kinematics and cross 
sections, 

• resonance decays that (more or less) form part of 
the hard subprocess (such as W, Z, t or h), 

• initial- and final-state QCD and QED show- 
ers (or, as an alternative, higher-order matrix 
elements, including a consistent treatment of 
virtual-correction terms) , 

• multiple parton-parton interactions, 

• beam remnants, 

• hadronization, 

• decay chains of unstable particles, and 

• general utility and analysis routines (such as jet 
finding) . 

However, even if a Monte Carlo includes all the physics 
we currently know of, there is no guarantee that 
not some important aspect of the physics is missing. 
Certain assumptions and phenomenological models 
inside the program are not well tested and will not 
necessarily hold when extrapolated to different energy 
regimes. For example, the strong-interaction dynamics 
in QCD remains unsolved and thereby unpredictable 
in an absolute sense. 

The PYTHIA 6.1 program was released in March 
1997, as a merger of JETSET 7.4, PYTHIA 5.7 
[ I and SPYTHIA [ ||. It addresses all of 
the aspects listed above. The current subversion is 
PYTHIA 6.136, which contains over 50,000 lines of 
Fortran 77 code. The code, manuals and sample main 
programs may be found at 

http : //www . thep . lu . se/'^torbj orn/Pythia . html . 

The two other programs of a similar scope are 
HERWIG[ I and IS A JET [ | 0. For parton-level 
processes, many more programs have been written. 

^http : //hepwww.rl . ac .uk/theory/seymour/herwig/ 
f tp : / /penguin . phy . bnl .gov/ pub/ isa j et 



The availability of several generators provides for use- 
ful cross-checks and a healthy competition. Since the 
physics of a complete hadronic event is very complex 
and only partially understood from first principles, 
one should not prematurely converge on one single 
approach. 

3.2. PYTHIA 6.1 Main Nevi^s 

Relative to previous versions, the main news in 
PYTHIA 6.1 includes 

• a renaming of the old JETSET program ele- 
ments to begin with PY, therefore now standard 
throughout, 

• new SUSY processes and improved SUSY simu- 
lation relative to SPYTHIA, and new PDG codes 
for sparticles, 

• new processes for Higgs (including doubly-char- 
ged in left-right symmetric models), technicolor, 
. . . , 

• several improved resonance decays, including an 
alternative Higgs mass shape, 

• some newer parton distributions, such as CTEQ5 

• initial-state showers matched to some matrix 
elements, 

• new options for final-state giuon splitting to a 
pair of c/b quarks and modified modeling of 
initial-state flavor excitation, 

• an energy-dependent p_Lmin in multiple interac- 
tions, 

• an improved modeling of the hadronization of 
small-mass strings, of importance especially for 
c/b, and 

• a built-in package for one-dimensional histo- 
grams (based on GBOOK). 

Some of these topics will be further studied below. 
Other improvements, of less relevance for pp colliders, 
include 

• improved modeling of gluon emission off c/b 
quarks in e^e^, 

• color rearrangement options for W~^W~ events, 

• a Bose-Einstein algorithm expanded with new 
options, 

• a new alternative baryon production scheme [ 

• QED radiation off an incoming muon, 

• a new machinery to handle real and virtual pho- 
ton fluxes, cross sections and parton distributions 
[1^, and 

• new standard interfaces for the matching to 
external generators of two, four and six fermions 
(and of two quarks plus two gluons) in e~^e~ . 
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The current list of over 200 different subprocesses 
covers topics such as hard and soft QCD, heavy flavors, 
DIS and 77, electroweak production of 7*/Z° and 
(singly or in pairs), production of a light or a heavy 
Standard Model Higgs, or of various Higgs states in 
supersymmetric (SUSY) or left-right symmetric mod- 
els, SUSY particle production (sfermions, gauginos, 
etc.), technicolor, new gauge bosons, compositeness, 
and leptoquarks. 

Needless to say, most users will still find that their 
particular area of interest is not as well addressed as 
could be wished. In some areas, progress will require 
new ideas, while lack of time and manpower is the 
limiting factor in others. 



3.3. Matching To Matrix Elements 

The matrix-element (ME) and parton-shower (PS) 
approaches to higher-order QCD corrections both have 
their advantages and disadvantages. The former offers 
a systematic expansion in orders of as, and a pow- 
erful machinery to handle multiparton configurations 
on the Born level, but loop calculations are tough 
and lead to messy cancellations at small resolution 
scales. Resummed matrix elements may circumvent 
the latter problem for specific quantities, but then do 
not provide exclusive accompanying events. Parton 
showers are based on an improved leading-log (almost 
next-to-leading-log) approximation, and so cannot be 
accurate for well separated partons, but they offer a 
simple, process-independent machinery that gives a 
smooth blending of event classes (by Sudakov form 
factors) and a natural match to hadronization. It is 
therefore natural to try to combine these descriptions, 
so that ME results are recovered for widely separated 
partons while the PS sets the subjet structure. 

For final-state showers in — > qq, corrections to 
the showering were considered quite a while ago [ Q , 
e.g. by letting the shower slightly overpopulate the qqg 
phase space and then using a Monte Carlo veto tech- 
nique to reduce down to the ME level. This approach 
easily carries over to showers in other color-singlet 
resonance decays, although the various relevant ME's 
have not all been implemented in PYTHIA so far. 

A similar technique is now available for the descrip- 
tion of initial-state radiation in the production of a 
single color-singlet resonance, such as j* /Z'^/W^ [ |39|] . 
The basic idea is to map the kinematics between the 
PS and ME descriptions, and to find a correction factor 
that can be applied to hard emissions in the shower 
so as to bring agreement with the matrix- element 
expression. Some simple algebra shows that, with 
the PYTHIA shower kinematics definitions, the two 



QQ 



gW^ emission rates disagree by a factor 



R. 



qq'^gW 



(do-/dt)ME 
(do-/dt) 



2to^s 



PS 



which is always between 1/2 and 1. The shower can 
therefore be improved in two ways, relative to the 
old description. Firstly, the maximum virtuality of 



emissions is raised from Q 



2 

max 



to Q 



2 

max 



S, 

i.e. the shower is allowed to populate the full phase 
space. Secondly, the emission rate for the final (which 
normally also is the hardest) q qg emission on each 
side is corrected by the factor R{s,t) above, so as 
to bring agreement with the matrix- element rate in 
the hard-emission region. In the backwards evolution 
shower algorithm [ this is the first branching 
considered. 

The other possible 0{as) graph is qg — > q'W^, 
where the corresponding correction factor is 



(d(T/dt)ME 
(do-/dt) 



PS 



(s — tti'Iy)'^ + rn' 



w 



which lies between 1 and 3. A probable reason for the 
lower shower rate here is that the shower does not ex- 
plicitly simulate the s-channel graph qg q* q'W. 
The g ^ qq branching therefore has to be preweighted 
by a factor of 3 in the shower, but otherwise the 
method works the same as above. Obviously, the 
shower will mix the two alternative branchings, and 
the correction factor for a final branching is based on 
the current type. 

The reweighting procedure prompts some other 
changes in the shower. In particular, u < trans- 
lates into a constraint on the phase space of allowed 
branchings. 

Our published comparisons with data on the W 
p±_ spectrum show quite a good agreement with this 
improved simulation [ |3^. A worry was that an 
unexpectedly large primordial k±, around 4 GeV, was 
required to match the data in the low-p_LVF region. 
However, at that time we had not realized that the data 
were not fully unsmeared. The required primordial fcj^ 
is therefore likely to drop by about a factor of two [ 

It should be noted that also other approaches to 
the same problem have been studied recently. The 
HERWIG one requires separate treatments in the 
hard- and soft-emission regions J_ 4l[ . Another, more 
advanced PYTHIA-based one [ ^|, also addresses the 
next-to-leading order corrections to the total W cross 
section, while the one outlined above is entirely based 
on the leading-order total cross section. There is also 
the possibility of an extension to Higgs production [ 
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Summarizing, we now start to believe we can handle 
initial- and final-state showers, with next-to-leading- 
order accuracy, in cases where these can be separated 
by the production of color singlet resonances — even if 
it should be realized that much work remains to cover 
the various possible cases. That still does not address 
the big class of QCD processes where the initial- 
and final-state radiation does not factorize. Possibly, 
correction factors to showers could be found also 
here. Alternatively, it may become necessary to start 
showers from given parton configurations of varying 
multiplicity and with virtual-correction weights, as 
obtained from higher-order ME calculations. So far, 
PYTHIA only implements a way to start from a given 
four-parton topology in e^e~ annihilation, picking one 
of the possible preceding shower histories as a way to 
set constraints for the subsequent shower evolution [ 
^ . This approach obviously needs to be extended in 
the future, to allow arbitrary parton configurations. 
Even more delicate will be the consistent treatment of 
virtual corrections [ where much work remains. 

3.4. Charm And Bottom Hadronization 

Significant asymmetries are observed between the 
production of D and D mesons in 7r~p collisions, with 
hadrons that share some of the ■n~ flavor content very 
much favored at large in the tt" fragmentation 
region [ . This behavior was qualitatively predicted 
by PYTHIA; in fact, the predictions were for somewhat 
larger effects than seen in the data. The new data has 
allowed us to go back and take a critical look at the 
uncertainties that riddle the heavy-flavor description [ 
Many effects are involved, and we limit ourselves 
here to mentioning only one. 

A hadronic event can be subdivided into sets of 
partons that form separate color singlets. These sets 
are represented by strings, that e.g. stretch from a 
quark end via a number of intermediate gluons to an 
antiquark end. The string has a mass, which can be 
calculated from the energy-momentum of its partons. 
Three different mass regions for the strings may be 
distinguished in the process of hadronization. 

1. Normal string fragmentation. This is the ideal 
situation, when each string has a large invariant 
mass, and the standard iterative fragmentation 
scheme [ Q works well. In practice, this 
approach can be used for all strings with a mass 
above a cut-off of a few GeV. 

2. Cluster decay. If a string is produced with a small 
invariant mass, then it is possible that only two- 
body final states are kinematically accessible. 
The traditional iterative Lund scheme is then 
not applicable. We call such a low-mass string 
a cluster, and treat it separately. In recent pro- 



gram versions, the modeling has been improved 
to give a smooth match onto the standard string 
scheme in the high-cluster-mass limit. 
3. Cluster collapse. This is the extreme case of 
the above situation, where the string mass is so 
small that the cluster cannot decay into even two 
hadrons. It is then assumed to collapse directly 
into a single hadron, which inherits the flavor 
content of the string endpoints. The original 
continuum of string/cluster masses is replaced 
by a discrete set of hadron masses. Energy 
and momentum then cannot be conserved inside 
the cluster, but must be exchanged with other 
objects within the local neighborhood. This 
description has also been improved. 
Because the mass of the charm and bottom par- 
tons are not negligible in the fragmentation process, 
the improved treatment of low-mass systems will 
have relatively more impact on charm and bottom 
hadronization. In general, flavor asymmetries are 
predicted to be smaller for bottom than for charm, 
and smaller at higher energies (except possibly at very 
large rapidities). Therefore, we do not expect any 
spectacular effects at the Tevatron. However, other 
nontrivial features of fragmentation may persist at 
higher energies, like a non-negligible systematic shift 
between the rapidity of a heavy quark parton and that 
of the hadron produced from it [ The possibility 
of such effects should be considered whenever trying 
to relate heavy flavor measurements to parton level 
calculations. 

3.5. Multiple Interactions 

Because of the composite nature of hadrons, several 
parton pairs may interact in a typical hadron-hadron 
collision [ |4^. Over the years, evidence for this 
mechanism has accumulated, such as the recent direct 
observation by CDF [ However, the occurrence of 
two hard interactions in one hadronic collision is just 
the tip of the iceberg. In the PYTHIA model, most 
interactions are at lower pj_ , where they are not visible 
as separate jets but only contribute to the underlying 
event structure. As such, they are at the origin of 
a number of key features, like the broad multiplicity 
distributions, the significant forward-backward multi- 
plicity correlations, and the pedestal effect under jets. 

Since the perturbative jet cross section is divergent 
for 0, it is necessary to regularize it, e.g. by a 

cut-off at some p_Lmin scale. That such a regularization 
should occur is clear from the fact that the incoming 
hadrons are color singlets — unlike the colored partons 
assumed in the divergent perturbative calculations — 
and that therefore the color charges should screen 
each other in the p± limit. Also other damping 
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mechanisms are possible [ Fits to data typically 
give p_Lmin ~ 2 GeV, which then should be interpreted 
as the inverse of some color screening length in the 
hadron. 

One key question is the energy-dependence of p_Linin; 
this may be relevant e.g. for comparisons of jet rates 
at different Tevatron energies, and even more for any 
extrapolation to LHC energies. The problem actually 
is more pressing now than at the time of our original 
study [ since nowadays parton distributions are 
known to be rising more steeply at small x than the flat 
xf{x) behavior normally assumed for small before 
HERA. This translates into a more dramatic energy 
dependence of the multiple-interactions rate for a fixed 

PJ_min . 

The larger number of partons also should increase 
the amount of screening, as confirmed by toy simula- 
tions [ |5^ . As a simple first approximation, p^min is 
assumed to increase in the same way as the total cross 
section, i.e. with some power e « 0.08 [ |5^ that, via 
reggeon phenomenology, should relate to the behavior 
of parton distributions at small x and Q^. Thus the 
new default in PYTHIA is 



(1.9 GcV) 



1 TeV 



0.08 



3.6. Interconnection Effects 

The widths of the W , Z and t are all of the order 
of 2 GeV. A Standard Model Higgs with a mass above 
200 GeV, as well as many supersymmetric and other 
"Beyond the Standard Model" particles would also 
have widths in the multi-GeV range. Not far from 
threshold, the typical decay times r^l/FwClfm^ 
"Tiad ~ 1 fm- Thus hadronic decay systems overlap, 
between a resonance and the underlying event, or 
between pairs of resonances, so that the final state may 
not contain independent resonance decays. 

So far, studies have mainly been performed in the 
context of W pair production at LEP2. Pragmatically, 
one may here distinguish three main eras for such 
interconnection: 

1. Perturbative: this is suppressed for gluon ener- 
gies w > F by propagator/timescale effects; thus 
only soft gluons may contribute appreciably. 

2. Non-perturbative in the hadroformation process: 
normally model-led by a color rearrangement be- 
tween the partons produced in the two resonance 
decays and in the subsequent parton showers. 

3. Non-perturbative in the purely hadronic phase: 
best exemplified by Bose-Einstein effects. 

The above topics are deeply related to the unsolved 
problems of strong interactions: confinement dyna- 
mics, 1/-/Vq effects, quantum mechanical interferences, 
etc. Thus they offer an opportunity to study the 



dynamics of unstable particles, and new ways to probe 
confinement dynamics in space and time [ |5j, 
but they also risk to limit or even spoil precision 
measurements . 

It is illustrative to consider the impact of inter- 
connection effects on the W mass measurements at 
LEP2. Perturbative effects are not likely to give 
any significant contribution to the systematic error, 
{5mw) ~ 5 MeV [ Color rearrangement is not 

understood from first principles, but many models 
have been proposed to model effects [ 57 , 

a conservative estimate gives {5mw) ^ 40 MeV 



and 
For 



Bose-Einstein again there is a wide spread in models, 
and an even wider one in results, with about the same 
potential systematic error as above [ 58, E^, Is^. The 



total QCD interconnection error is thus below TOjt 
in absolute terms and 0.1% in relative ones, a small 
number that becomes of interest only because we aim 
for high accuracy. 



A study of e+e- ^ tt ^ bW+bW- 



bbe+vd'^v' 



near threshold gave a realistic interconnection uncer- 
tainty of the top mass of around 30 MeV, but also 
showed that slight mistreatment of the combined color 
and showering structure could blow up this error by 
a factor of ten [ For hadronic top decays, errors 
could be much larger. 

The above numbers, when applied to hadronic 
physics, are maybe not big enough to cause an im- 
mediate alarm. The addition of a colored underlying 
event — with a poorly-understood multiple-interaction 
structure as outlined above — has not at all been 
considered so far, however, and can only make matters 
worse in hadronic physics than in e+e^. This is clearly 
a topic for the future, where we should be appropri- 
ately humble about our current understanding, at least 
when it comes to performing precision measurements. 

QCD interconnection may also be at the root of 
a number of other, more spectacular effects, such as 
rapidity gaps and the whole Pomeron concept [|6|] , and 
the unexpectedly large rate of quarkonium production 

[M- 

3.7. The Future: On To C++ 

Finally, a word about the future. PYTHIA continues 
to be developed. On the physics side, there is a 
need to increase the support given to different physics 
scenarios, new and old, and many areas of the general 
QCD machinery for parton showers, underlying events 
and hadronization require further improvements, as we 
have seen. 

On the technical side, the main challenge is a 
transition from Fortran to C++, the language of choice 
for Run II (and LHC). To address this, the PYTHIA 7 
project was started in January 1998, with L. Lonnblad 
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bearing the main responsibility. A similar project, but 
more ambitious and better funded, is now starting 
up for HERWIG, with two dedicated postdoc-level 
positions and a three-year time frame. 

For PYTHIA, what exists today is a strategy doc- 
ument [ and code for the event record, the 
particle object, some particle data and other data base 
handling, and the event generation handler structure. 
All of this is completely new relative to the Fortran 
version, and is intended to allow for a much more 
general and flexible formulation of the event gener- 
ation process. The first piece of physics, the string 
fragmentation scheme, is being implemented by M. 
Bertini, and is nearing completion. The subprocess 
generation method is being worked on for the simple 
case of e+e~ — > ^ qq. The hope is to have a "proof 
of concept" version soon, and some of the current 
PYTHIA functionality up and running by the end of 
2000. It will, however, take much further effort after 
that to provide a program that is both more and better 
than the current PYTHIA 6 version. It is therefore 
unclear whether PYTHIA 7 will be of much use during 
Run II, except as a valuable exercise for the future. 

4. A Comparison of the Predictions from 
Monte Carlo Programs and Transverse Mo- 
mentum Resummation 

by C. Baldzs, J. Huston, I. Puljak, S. Mrenna 

4.1. Introduction 

Monte Carlo programs including parton showering, 
such as PYTHIA[|, HERWIG[|l| and ISAJET[|, 
are commonly used by experimentalists, both as a way 
of comparing experimental data to theoretical predic- 
tions, and also as a means of simulating experimental 
signatures in kinematic regimes for which there is no 
experimental data (such as that appropriate to the 
LHC). The final output of the Monte Carlo programs 
consists of the 4- vectors of a set of stable particles (e.g., 
e, /i, TT, 7); this output can either be compared to recon- 
structed experimental quantities or, when coupled with 
a simulation of a detector response, can be directly 
compared to raw data taken by the experiment, and/or 
passed through the same reconstruction procedures as 
the raw data. In this way, the parton shower programs 
can be more useful to experimentalists than analytic 
calculations performed at high orders in perturbation 
theory. Indeed, almost all of the physics plots in the 
ATLAS physics TDR [ lOS] involve comparisons to 



PYTHIA(version 5.7). 

Here, we are concerned with the predictions of 
parton shower Monte Carlo programs and those from 
certain analytic calculations which resum logarithms 
associated with the transverse momentum of partons 



initiating the hard scattering. Most analytic cal- 
culations of this kind are either based on or origi- 
nate from the formalism developed by J. Collins, D. 
Soper, and G. Sterman (CSS), which we choose as 
the analytic "benchmark" of this section. Both the 
parton showering and analytic calculations describe 
the effects of multiple soft gluon emission from the 
incoming partons, which can have a profound effect 
on the kinematics of gauge or Higgs bosons and their 
decay products produced in hadronic collisions. This 
may have an impact on the signatures of physics 
processes at both the trigger and analysis levels, and 
thus it is important to understand the reliability of 
such predictions. The best method for testing the 
reliability is a direct comparison of the predictions 
to experimental data. If no experimental data is 
available, then some understanding of the reliability 
may by gained by simply comparing the predictions of 
different calculational methods. 

4.2. Parton Showering and Resummation 

Parton showering is the backwards evolution of an 
initial hard scattering process, involving only a few 
partons at a high scale Q^^^x reflecting large virtuality, 
into a complicated, multi-parton configuration at a 
much lower scale Q^j^ typical of hadronic binding ener- 
gies. In practice, one does not calculate the probability 
of arriving at a specific multi-parton configuration 
all at once. Instead, the full shower is constructed 
in steps, with evolution down in virtuality with 
no parton emission, followed by parton emission, and 
then a further evolution downward with no emission, 
etc., until the scale Q^j^ is reached. The essential 
ingredient for this algorithm is the probability of 
evolving down in scale with no parton emission or 
at least no resolvable parton emission. This can be 
derived from the DGLAP equation for the evolution 
of parton distribution functions. One finds that the 
probability of no emission P equals 1 — exp(— 5), where 
S is the Sudakov form factor, a function of virtuality 
and the momentum fraction x carried by a parton. 

A key ingredient in the parton showering algorithm 
is the conservation of energy-momentum at every step 
in the cascade. The transverse momentum of the 
final system partly depends on the opening angle 
between the mother and daughter partons in each 
emission. Furthermore, after each emission, the entire 
multi-parton system is boosted to the center-of-mass 
frame of the two virtual partons, until at the end of the 
shower one is left with two primordial partons which 
are on the mass shell and essentially parallel with the 
incoming hadrons. These boosts also influence the final 
transverse momentum. 

Parton showering resums primarily the leading log- 
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arithms - those resummed by the DGLAP equations 
- which are universal, i.e. process independent, and 
depend only on the given initial state. In this lies 
one of the strengths of the parton shower approach, 
since it can be incorporated into a wide variety of 
physical processes. An analytic calculation, in compar- 
ison, can resum many other types of potentially large 
logarithms, including process dependent ones. For 
example, the CSS formalism in principle sums all of the 
logarithms with lv\ in their arguments, where, for 
the example of Higgs boson production, Q is the four 
momentum of the Higgs boson and pt is its transverse 
momentum. All of the "dangerous logs" are included in 
the Sudakov exponent, which can be written in impact 
parameter (6) space as: 



2 

S{Q,b)=l % 



A {as{jl)) In ( ^ ) +B {as[jl)) 



with the A and B functions being free of large 
logarithms and calculable in fixed-order perturbation 
theory: 



oo 



(17) 



These functions contain an infinite number of coeffi- 
cients, with the A^"^ being universal to a given initial 
state, while the i?^") are process dependent. In prac- 
tice, the number of towers of logarithms included in the 
Sudakov exponent depends on the level to which a fixed 
order calculation was performed for a given process. 
For example, if only a next-to- leading order calculation 
is available, only the coefficients A^^^ and B'^^^ can 
be included. If a NNLO calculation is available, then 
and ^(2) can be extracted and incorporated into 
a resummation calculation, and so on. This is the case, 
for example, for boson production. So far, only the 
^(2) and coefficients are known for Higgs 
production, but the calculation of S^^^ is in progress [ 
|109| . If we try to interpret parton showering in the 
same language, then we can say that the parton shower 
Sudakov exponent always contains a term analogous to 
A'^^\ It was shown in Reference [ |110| that a suitable 
modification of the Altarelli-Parisi splitting function, 
or equivalently the strong coupling constant as, also 
effectively approximates the A^'^^ coefficient .|^ 

In contrast with parton showering, analytic resum- 
mation calculations integrate over the kinematics of 
the soft gluon emission, with the result that they are 



limited in their predictive power. While the parton 
shower maintains an exact treatment of the branching 
kinematics, the original CSS formalism imposes no 
kinematic penalty for the emission of the soft gluons, 
although an approximate treatment of this can be 
incorporated into a numerical implementation, like 
ResBos 
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Neither parton showering nor ana- 
lytic resummation reproduces kinematic configurations 
where one hard parton is emitted at large pt- In 
the parton shower, matrix element corrections can be 
imposed [ while, in the analytic resummation 

calculation, matching is necessary. 

With the appropriate input from higher order cross 
sections, a resummation calculation has the corre- 
sponding higher order normalization and scale depen- 
dence. The normalization and scale dependence for 
the Monte Carlo, though, remains that of a leading 
order calculation - though see Ref. [ ^ and the related 
contribution to these proceedings for an idea of how to 
include these at NLO. The parton showering occurs 
with unit probability after the hard scattering, so it 
does not change the total cross section.^ 

Given the above discussion, one quantity which 
should be well-described by both calculations is the 
shape of the transverse momentum [pr) distribution of 
the final state electroweak boson in a subprocess such 
as qq WX, ZX or gg — > HX, where most of the 
Pt is provided by initial state parton showering. The 
parton showering supplies the same sort of transverse 
kick as the soft gluon radiation in a resummation 
calculation. Indeed, very similar Sudakov form factors 
appear in both approaches, with the caveats about the 
A^"^ and i?'-"-' terms mentioned previously. 

At a point in its evolution corresponding to a 
virtuality on the order of a few GeV, the parton shower 
is stopped and the effects of gluon emission at softer 
scales must be parameterized and inserted by hand. 
Typically, a Gaussian probability distribution function 
is used to assign an extra "primordial" kr to the 
primordial partons of the shower (the ones which are 
put on the mass shell at the end of the backwards 
showering). In PYTHIA, the default is a constant 
value of kr- Similarly, there is a somewhat arbitrary 
division between perturbative and non-perturbative 
regions in a resummation calculation. Sometimes 
the non-perturbative effects are also parametrized by 
Gaussian distributions in b or Qt space. In general, the 
value for the non-perturbative (fcr) needed in a Monte 
Carlo program will depend on the particular kinemat- 



^This is rigorously true only for the high parton x or y/r region. 



^Technically, one could add the branching for q — > q+Higgs in 
the shower, which would have the capability of increasing some- 
what the Higgs cross section; however, the main contribution to 
the higher order /sT-factor comes from the virtual corrections and 
the 'Higgs Bremsstrahlung' contribution is negligible. 
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ics being investigated. In the case of the resummation 
calculation the non-perturbative physics is determined 
from fits to fixed target data and then automatically 
evolved to the kinematic regime of interest. 

A value for the average non-perturbative kx of 
greater than 1 GeV does not imply that there is an 
anomalous intrinsic kx associated with the parton 
size; rather this amount of (kx) needs to be supplied 
to provide what is missing in the truncated parton 
shower. If the shower is cut off at a higher virtuality, 
more of the "non-perturbative" kx will be needed. 

4.3. Boson Production at the Tevatron 

The 4-vector of a boson, and thus its transverse 
momentum, can be measured with great precision in 



the 



decay mode. Resolution effects are relatively 



minor and are easily corrected for. Thus, the Z'^ 
px distribution is a great testing ground for both 
the resummation and Monte Carlo formalisms for soft 
gluon emission. The corrected px distribution for Z'^ 
bosons in the low px region for the CDF experimenij^ 
is shown in Figure ^, compared to both the resummed 
prediction from ResBos, and to two predictions from 
PYTHIA (version 6.125). One PYTHIA prediction uses 
the default (rms)|^ value of intrinsic kx of 0.44 GeV and 
the second a value of 2.15 GeV per incoming parton.^ 
The latter value was found to give the best agreement 
between PYTHIA and the dataj^ AU of the predictions 
use the CTEQ4M parton distributions [|ll2). The shift 
between the two PYTHIA predictions at low px is 
clearly evident. As might have been expected, the high 
Px region (above 10 GeV) is unaffected by the value 
of the non-perturbative kx- Note the kx imparted to 
the incoming partons at their lowest virtuality, Qq, is 
greatly reduced in its effect on the Z'^ px distribution. 
This dilution arises because the center-of-mass energy 
of the "primordial" partons is typically much larger 
than that of the original hard scattering. Therefore, 
the transverse /3 of the boost applied to the Z'^ boson 
to transform it to the frame where the "primordial" 
partons have transverse momentum kx is small. 

As an exercise, one can transform the resummation 
formula in order to bring it to a form where the non- 
perturbative function acts as a Gaussian type smearing 
term. Using the Ladinsky-Yuan parameterization [ 
L14| of the non-perturbative function in ResBos leads 



'^We thank Willis Sakumoto for providing the figures for 

production as measured by CDF 

*For a Gaussian distribution fc™= = 1.13(fcT>. 

previous publication [ b9| indicated the need for a substan- 
tially larger non-perturbative (kx), of the order of 4 GeV for 
the case of W production at the Tevatron. The data used in the 
comparison, however, were not corrected for resolution smearing, 
a fairly large effect for the case ofW^eu production and decay. 
^"A similar conclusion has been re ached for comparisons of the 
CDF Z° pT data with HERWIG. 



to an rms value for the effective kx smearing param- 
eter, for production at the Tevatron, of 2.5 GeV. 
This is similar to that needed for PYTHIA and HERWIG 
to describe the Z° production data at the Tevatron. 

In Figure ^, the normalization of the resummed 
prediction has been rescaled upwards by 8.4%. The 
PYTHIA prediction was rescaled by a factor of 1.3-1.4 
(remember that this is only a leading order compari- 
son) for the shape comparison. 
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Figure 2. The Z^ px distribution (at low px) from 
CDF for Run 1 compared to predictions from ResBos 
and from PYTHIA. The two PYTHIA predictions use the 
default (rms) value for the non-perturbative kx (0.44 
GeV) and the value that gives the best agreement with 
the shape of the data (2.15 GeV). 



As stated previously, the resummed prediction cor- 
rectly describes the shape of the Z^ px distribution at 
low Px, although there is still a noticeable difference 
in shape between the Monte Carlo and the resummed 
prediction. It is interesting to note that if the process 
dependent coefficients {B^^^ and B^"^^) were not incor- 
porated into the resummation prediction, the result 
would be an increase in the height of the peak and a 
decrease in the rate between 10 and 20 GeV, leading to 
a better agreement with the PYTHIA prediction 
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The PYTHIA and ResBos predictions both describe 
the data well over a wider px range than shown in the 
figure. Note especially the agreement of PYTHIA with 
the data at high px, made possible by explicit matrix 
element corrections (from the subprocesses qq —^ Z^g 
and gq — > Z^q) to the Z^ production process.^ 

Slightly different techniques are used for the matrix element 
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4.4. Diphoton Production 

Most of the comparisons between resummation cal- 
culations/Monte Carlos and data have been performed 
for Drell-Yan production, i.e. qq initial states. It is 
also interesting to examine diphoton production at the 
Tevatron, where a large fraction of the contribution 
at low diphoton mass is due to gg scattering. The 
prediction for the di-photon kx distribution at the 
Tevatron, from PYTHIA (version 6.122), is shown in 
Figure ^, using the experimental cuts applied in the 



CDF analysis [ |ll6 | 



Tevatron, = 1.8 TcV, PYTHIA 6.122 
lri^<l,PT,>12GeV 

gg + qq — > yy + X, a = 1 2.08 pb 

gg-> 7y+ X, o = 4.90 pb 

qq— >YY + X,a= 7.18 pb 
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Figure 3. A comparison of the PYTHIA predictions 
for di-photon production at the Tevatron for the two 
different subprocesses, qq and gg. The same cuts are 
applied to PYTHIA as in the CDF di-photon analysis. 



corrections by PYTHIA [ ^ and by HERWIG [ |4l| . In PYTHIA, 
the parton shower probability distribution is applied over the 
whole phase space and the exact matrix element corrections are 
applied only to the branching closest to the hard scatter. In 
HERWIG, the corrections are generated separately for the regions 
of phase space unpopulated by HERWIG (the 'dead zone') and the 
populated region. In the dead zone, the radiation is generated 
according to a distribution using the first order matrix element 
calculation, while the algorithm for the already populated region 
applies matrix element corrections whenever a branching is 
capable of being 'the hardest so far'. 



It is interesting to note that about half of the 
di-photon cross section at the Tevatron is due to the 
gg subprocess, and that the di-photon px distribution 
is noticeably broader for the gg subprocess than the qq 
subprocess. The gg subprocess predictions in ResBos 
agree well with those from PYTHIA while the qq pT 
distribution is noticeably broader in ResBos. The 
latter behavior is due to the presence of the Y piece 
(fixed-order corrections) in ResBos at moderate px, 
and the matching of the qq cross section to the fixed 
order qq — > 773 at high pT- The corresponding matrix 
element correction is not in PYTHIA. It is interesting 
to note that the PYTHIA and ResBos predictions for 
gg 77 agree in the moderate pT region, even though 
the ResBos prediction has the Y piece present and 
is matched to the matrix element piece gg 77(7 
at high pt, while there is no such matrix element 
correction for PYTHIA. This shows that the Y piece 
correction is not important for the gg subprocess, 
which is the same conclusion that was reached in Ref. \ 



117] . This is probably a result of steep decHne in the gg 



parton-parton with increasing partonic center of mass 
energy, v^. This falloff tends to suppress the size of the 
Y piece since the production of the di-photon pair at 
higher pr requires larger xi, X2 values. In the default 
CSS formalism, there is no such kinematic penalty 
in the resummed piece since the soft gluon radiation 
comes for "free." (Larger xi and X2 values are not 
required.) 

A comparison of the CDF di-photon data to NLO 
[ |11^ and resummed (ResBos) QCD predictions has 
been performed, but the analysis is still in progress, so 
the results are not presented here. The transverse mo- 
mentum distribution, in particular, is sensitive to the 
effects of the soft gluon radiation and better agreement 
can be observed with the ResBos prediction than with 
the NLO one. A much more precise comparison with 
the effects of soft gluon radiation will be possible with 
the 2 fb~^ or greater data sample that is expected for 
both CDF and D0 in Run 2. 

4.5. Higgs Boson Production 

A comparison of the two versions of PYTHIA and 
of ResBos is shown in Figure ^ for the case of the 
production of a Higgs boson with mass 100 GeV at the 
Tevatron with center-of-mass energy of 2.0 TeV. The 
same qualitative features are observed at the LHC: the 
newer version of PYTHIA agrees better with ResBos in 
describing the low pT shape, and there is a falloff at 
high pt unless the hard scale for showering is increased. 
The default (rms) value of the non-perturbative kx 
(0.44 GeV) was used for the PYTHIA predictions. Note 
that the peak of the resummed distribution has moved 
to Pt ~ 7 GeV (compared to about 3 GeV for Z° 
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production at the Tevatron). This is due primarily 
to the larger color factors associated with initial state 
gluons {Ca — 3) rather than quarks {Cp = 4/3). 

The newer version of PYTHIA agrees well with ResBos 
at low to moderate pr, but falls below the resummed 
prediction at high pT- This is easily understood: 
ResBos switches to the NLO Higgs + jet matrix 
element at high pT while the default PYTHIA can 
generate the Higgs pr distribution only by initial 
state gluon radiation, using as default a maximum 
scale equal to the Higgs boson mass. High pT Higgs 
boson production is another example where a, 2 —f 1 
Monte Carlo calculation with parton showering can 
not completely reproduce the exact matrix element 
calculation without the use of matrix element correc- 
tions. The high pT region is better reproduced if the 
maximum virtuality Q^ax is set equal to the collider 
center-of-mass energy, s, rather than subprocess s. 
This is equivalent to applying the parton shower to 
all of phase space. However, the consequence is that 
the low pt region is now depleted of events, since the 
parton showering does not change the total production 
cross section. The appropriate scale to use in PYTHIA 
(or any Monte Carlo) depends on the pT range to 
be probed. If matrix element information is used 
to constrain the behavior, the correct high pT cross 
section can be obtained while still using the lower 
scale for showering. The incorporation of matrix 
element corrections to Higgs production (involving the 
processes gq qH,qq gH, gg gH) is the next 
logical project for the Monte Carlo experts, in order to 
accurately describe the high pT region. 

The older version of PYTHIA produces too many 
Higgs events at moderate pr (in comparison to Res- 
Bos) at both the Tevatron and the LHC. Two changes 
have been implemented in the newer version. The first 
change is that a cut is placed on the combination of z 
and values in a branching: u = / z — s{l — z) < 0, 
where s refers to the subsystem of the hard scattering 
plus the shower partons considered to that point. 
The association with u is relevant if the branching is 
interpreted in terms of a 2 ^ 2 hard scattering. This 
requirement is not fulfilled when the value of the 
space-like emitting parton is little changed and the z 
value of the branching is close to unity. This affects 
mainly the hardest emission (largest Q'^). The net 
result of this requirement is a substantial reduction 
in the total amount of gluon radiation [ |119| . Such 
branchings are kinematically allowed, but since matrix 
element corrections would assume initial state partons 
to have = 0, a non-physical u results (and thus 
no possibility to impose matrix element corrections). 
The correct behavior is beyond the predictive power of 
leading log Monte Carlos. 



^ 0.04 
.S 0.035 
^ 0.03 
, 0.025 



■§ 0.015 
0.01 
0.005 




•ft 




gg — > H + X at Tevatron 

m„ = 100 GeV, CTEQ4M, Vs = 2 TeV 

ResBos 98.07.14 

PYTHIA 5.7 default 

in THLA 122 defaull 

PYTHIA 6.122 q'„ 



40 45 50 
Pj.(GeV) 



D 
"a 




80 90 100 

Pj.(GeV) 



Figure 4. A comparison of predictions for the Higgs 
PT distribution at the Tevatron from ResBos and from 
two recent versions of PYTHIA. The ResBos and PYTHIA 
predictions have been normalized to the same area. 



In the second change, the parameter for the min- 
imum gluon energy emitted in space-like showers is 
modified by an extra factor roughly corresponding to 
the 1/7 factor for the boost to the hard subprocess 
frame [ 119|. The effect of this change is to increase 
the amount of gluon radiation. Thus, the two effects 
are in opposite directions but with the first effect being 
dominant. 

This difference in the px distribution between the 
two versions of PYTHIA could have an impact on the 
analysis strategies for Higgs searches at the LHC. 
For example, for the CMS detector, the higher pT 
activity associated with Higgs production in version 5.7 
would have allowed for a more precise determination 
of the event vertex from which the Higgs (decaying 
into two photons) originated. Vertex pointing with the 
photons is not possible in CMS, and the large number 
of interactions occurring with high intensity running 
will mean a substantial probability that at least one of 
the interactions will produce jets at low to moderate 
Et- [ 120 In principle, this problem could affect the 
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Pt distribution for all PYTHIA processes. In practice, 
the effect has manifested itself only in gg initial states, 
due to the enhanced branching probability. 

As an exercise, an 80 GeV W and an 80 GeV 
Higgs were generated at the Tevatron using PYTHIA5.7 [ 
[L2l| . A comparison of the distribution of values of u 
and the virtuality Q for the two processes indicates a 
greater tendency for the Higgs virtuality to be near the 
maximum value and for there to be a larger number of 
Higgs events with positive u (than W events). 

4.6. Comparison with HERWIG 

The variation between versions 5.7 and 6.1 of PYTHIA 
gives an indication of the uncertainties due to the types 
of choices that can be made in Monte Carlos. The 
requirement that u be negative for all branchings is a 
choice rather than an absolute requirement. Perhaps 
the better agreement of version 6.1 with ResBos is 
an indication that the adoption of the u restrictions 
was correct. Of course, there may be other changes 
to PYTHIA which would also lead to better agreement 
with ResBos for this variable. 

Since there are a variety of choices that can be made 
in Monte Carlo implementations, it is instructive to 
compare the predictions for the px distribution for 
Higgs boson production from ResBos and PYTHIA with 
that from HERWIG (version 5.6, also using the CTEQ4M 
parton distribution functions). The HERWIG prediction 
is shown in Figure || along with the PYTHIA and ResBos 
predictions, all normalized to the ResBos prediction. 
In all cases, the CTEQ4M parton distribution was 
used. The predictions from HERWIG and PYTHIA 6.1 
are very similar, with the HERWIG prediction matching 
the ResBos shape somewhat better at low px- 

4.7. Non-perturbative fc^ 

A question still remains as to the appropriate value 
of non-perturbative kx to input in the Monte Carlos 
to achieve a better agreement in shape, both at the 
Tevatron and at the LHC. Figure ^ compares the 
ResBos and PYTHIA predictions for the Higgs boson px 
distribution at the Tevatron. The PYTHIA prediction 
(now version 6.1 alone) is shown with several values 
of non-perturbative kx- Surprisingly, no difference is 
observed between the predictions with the different 
values of kx, with the peak in PYTHIA always being 
somewhat below that of ResBos. This insensitivity can 
be understood from the plots at the bottom of the two 
figures which show the sum of the non-perturbative 
initial state kx {kxi+kx2) at Qo and at the hard 
scatter scale Q. Most of the kx is radiated away, 
with this effect being larger (as expected) at the 

^■^The normalization factors (RcsBos/Monte Carlo) are PYTHIA 
(both versions){1.61) and HERWIG (1.76). 
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Figure 5. A comparison of predictions for the Higgs 
px distribution at the LHC from ResBos, two recent 
versions of PYTHIA and HERWIG. The ResBos, PYTHIA 
and HERWIG predictions have been normalized to the 
same area. 



LHC. The large gluon radiation probability from a 
gluon-gluon initial state (and the greater phase space 
available at the LHC) lead to a stronger degradation 
of the non-perturbative kx than was observed with 
production at the Tevatron. 

4.8. Conclusions 

An understanding of the signature for Higgs boson 
production at either the Tevatron or LHC depends 
upon the understanding of the details of soft gluon 
emission from the initial state partons. This soft gluon 
emission can be modeled either in a Monte Carlo or 
in a kx resummation program, with various choices 
possible in both implementations. A comparison of the 
two approaches is useful to understand the strengths 
and weaknesses of each. The data from the Tevatron 
that either exists now, or will exist in Run 2, will be 
extremely useful to test both approaches. 
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Figure 6. (top) A comparison of the PYTHIA predic- 
tions for the pT distribution of a 100 GeV Higgs at the 
Tevatron using the defauh (rnis) non-perturbative kx 
(0.44 GeV) and a larger value (4 GeV), at the initial 
scale Qo and at the hard scatter scale. Also shown is 
the ResBos prediction (bottom) The vector sum of the 
intrinsic kx (fcTi+fcT2) for the two initial state partons 
at the initial scale Qo and at the hard scattering scale 
for the two values of intrinsic kx ■ 
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5. MCFM: a parton-level Monte Carlo at NLO 
Accuracy 

by John Campbell and R.K. Ellis 

5.1. Introduction 

In Run II, experiments at the Tevatron will be 
sensitive to processes occurring at the femtobarn level. 
Of particular interest are processes which involve heavy 
quarks, leptons and missing energy, since so many of 
the signatures for physics beyond the standard model 
produce events containing these features. We have 
therefore written the program MCFM [ 123, |7^ which 
calculates the rates for a number of standard model 



processes. These processes are included beyond the 
leading order in the strong coupling constant where 
possible; in QCD this is the first order in which the nor- 
malization of the cross sections is determined. Because 
the program produces weighted Monte Carlo events, 
we can implement experimental cuts allowing realistic 
estimates of event numbers for an ideal detector con- 
figuration. MCFM is expected to give more reliable 
results than parton shower Monte Carlo programs, 
especially in phase space regions with well separated 
jets. On the other hand it gives little information 
about the phase space regions which are dominated 
by multiple parton emission. In addition, because the 
final state contains partons rather than hadrons, a full 
detector simulation cannot be performed directly using 
the output of MCFM. 

The processes already included in MCFM at NLO 
are as follows {Hi,H2 = p or p), 
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The decays of vector bosons and/or Higgs bosons are 
included. We have also included the leptonic decays of 
the T-lepton. As described below the implementation 
of NLO corrections requires the calculation of both the 
amplitude for real radiation and the virtual corrections 
to the Born level process. We have extensively used 
the one loop results of Bern, Dixon, Kosower et al. [ 
124], [ |125[ to obtain the virtual corrections to above 



processes. 

A future development path for the program would 
be to include the following processes at NLO: 



• Hi 

• Hi 



H2 
H2 



+ 2 jets 
Z + 2 jets . 
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In addition there are an number of processes which 
we have included only at leading order. This restriction 
to leading order is both a matter of expediency and be- 
cause the theoretical framework for including radiative 
corrections to processes involving massive particles is 



not yet complete 
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H, Z and top quark decays are included. 

5.2. General structure 

In order to evaluate the strong radiative corrections 
to a given process, we have to consider Feynman 
diagrams describing real radiation, as well as the 
diagrams involving virtual corrections to the tree level 
graphs. The corrections due to real radiation are dealt 
with using a subtraction algorithm[ 126 as formulated 
by Catani and Seymour [ 127|. This algorithm is 



based on the fact that the singular parts of the QCD 
matrix elements for real emission can be singled out 
in a process-independent manner. By exploiting this 
observation, one can construct a set of counter-terms 
that cancel all non-integrable singularities appearing 
in real matrix elements. The NLO phase space 
integration can then be performed numerically in four 
dimensions. 

The counter-terms that were subtracted from the 
real matrix elements have to be added back and 
integrated analytically over the phase space of the 
extra emitted parton in n dimensions, leading to poles 
in e = (n — 4)/2. After combining those poles with the 
ones coming from the virtual graphs, all divergences 
cancel, so that one can safely perform the limit e —^ 
and carry out the remaining phase space integration 
numerically. 

As an example of this procedure we consider the 
production of an on-shell W boson decaying to a 
lepton-antilepton pair. 



q{pi) + q{p2) ^ w+{u{p3) + e+(p4)), 

P2=PZ+ PA, {P3 + PaY = Ml 



(18) 



Pi 

In this case 

calculable analytically in 0{as) 
shows the result calculated in the MS scheme. The 
virtual corrections to (18) are of the Drell-Yan type 



the W boson rapidity distribution is 

Fig. I 
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Figure 7. The rapidity distribution for production 
in pp collisions at ^/s = 2TeV. 



and are well known [ 128|. They are expressible as 
an overall factor multiplying the lowest order matrix 
element squared, 



27r 



1 



r(l-e) 
2 3 



6 + 7r" 



(19) 



and must be combined with the real radiation contribu- 
tion. For example, gluon radiation from the qq initial 
state yields the subprocess 

qipi) + q{P2) Wivips) -h e+{pi)) + 9{Pb), 

Pl+P2^P3+P4+P5- (20) 

To eliminate the singular part of this subprocess, we 
generate a counter event with the kinematics of the 
2^2 process as follows 

q{xaPi) + q{P2) W{v{p3) + e+{pi)), 

XaPl +P2 =P3+P4: (21) 

where a Lorentz transformation has been performed on 
all j final state momenta 



Pi 



(22) 



such that pj — > pj for p^ collinear or soft. Thus 
the energy of the emitted gluon p^ is absorbed by pi , 
and the momentum components are absorbed by the 
transformation of the final state vectors. The phase 
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space has a convolution structure, 



dxd^'^'^\p4,P3;p2,Pi) X [dp5{pi,p2,x)] 

(23) 



where 



[dp5{pi,p2,: 



d^p^ 



j^^-^^5+{pi)Q{x)Q{l~x)5{x-Xa) (24) 

This phase space may be used to integrate out the 
dipole term D^^'^ , which is chosen to reproduce the 
singularities in the real matrix elements as the giuon 
(5) becomes soft or collinear to the quark (1), 
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Pi- Pb \l- Xa 

Performing the integration yields, 

/ dxD^'''^[dp5{puP2.x)] = 
Jo 

asCp ( X 
2ti \2pi-p2j r(l-e) 

-Pqq{x) + '^(1 - a:) I ^ + — - — 



(25) 



+2(1 + x2) 



log(l - x) 

l~x 



(26) 



with the Altarelli-Parisi function Pqq{x) given by 



Pqqix) 



{1-xy 



l-x+-S{l-x) (27) 



In order to obtain the complete counter-term, one must 
add the (identical) contribution from the dipole con- 
figuration D^^'^ that accounts for the gluon becoming 
collinear with the anti-quark. In a more complicated 
process, we would sum over a larger number of distinct 
dipole terms involving partons both in the initial and 
final states. In this simple case, we find the total 
counter-term contribution to the qq cross-section to be 



^CT 



UsCp f 47r/i 



2\ « 



2tt \ J r(l - e) 



2 .X N / 2 3 TT^ 

--Pqqix) +S{1-X) 1 ^ + --y 



-2pqq{x) logX + 4(l +x2) 



log(l - x) 



1-X 



where each of these terms leads to a different type of 
contribution in MCFM. The first term, proportional to 
Pqq{x), is canceled by mass factorization, up to some 
additional finite {0{e°)) pieces. The terms multiplying 
the delta- function S{1 — x) manifestly cancel the poles 
generated by the virtual graphs, given in equation (19), 
leaving an additional tt^ contribution. The remaining 
terms, which don't have the structure of the virtual 
contribution, are collected together and added sepa- 
rately in MCFM. 




Figure 8. The ratio of the contributions to the rapidity 
distribution of production 



In Fig. ^ we have plotted the three contributions 
to the W rapidity calculated using MCFM. The 
three contributions are a) the contribution of (real- 
counterterni) [the lower curve], b) the contribution 
of leading order + virtual -I- integrated counter-term 
[the upper-most curve] and c) the total contribution. 
All three terms have been normalized to the 0{as) 
rapidity distribution shown in Fig. |^. We see that 
(&), the leading order term, combined with the virtual 
correction and the results from the counterterm pro- 
vides the largest contribution to the cross section. The 
total contribution is a horizontal line at unity, showing 
the agreement between MCFM and the analytically 
calculated result. Only at the boundaries of the 
phase space at large y can the contribution of the real 
emission minus the counterterm become sizeable. 

5.3. Examples of MCFM results 

We first detail the input parameters used in our 
phenomenological estimates. The electroweak theory is 
specified by four numbers, MwTMz,a{Mz), and Gf, 
the values of which are given in Table ||, together with 
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Table 1 

Input parameters 



Mz,rz 


91.187,2.49 GeV 


Mw, 


80.41, 2.06 GeV 


mt.Tt 


175, 1.4 GeV 


a{Mz) 


1/128.89 


Gp 


1.16639 X 10"^ 


• 2 n 

sm uw 


0.228534483 


T r 

Viid 


U. 97500 




0.22220 


Vcd 


0.22220 


Yes 


0.97500 


Higgs mass (GeV) 


BR(iJ ^ 66) 


100 


0.8119 


110 


0.7697 


120 


0.6778 


130 


0.5254 



other necessary constants. Other derived parameters 
are e^gw and sin^ which, when defined as below, 
are effective parameters including the leading effects 
of top quark loops [ |130|. We use the the first 



of the MRS99 parton distributions [ 131| which has 
as[Mz) = 0.1175. 



9w 



= 47ra(Afz) 
= 4\/2GfM^ 

^2 



W-Bw = ^ (28) 
Tabic H shows the production cross sections for 



Table 2 

Diboson cross sections (in pb) at the Tevatron and the 
LHC 







a{W+Z) 1 a{W-Z) 


a{ZZ) 


2 TeV (pp) 


12.2 


2.02 


1.75 


14 TeV (pp) 


103.6 


27.2 1 17.7 


16.7 



di-boson production calculated using MCFM for 
collisions at -y/i — 2 TeV and for collisions at 
•^s — 14 TeV. The next-to-leading order corrections 
vary between approximately 30% and 50% of leading 
order and are almost entirely due to the virtual graphs. 
The numbers here are slightly different than the results 
in [ ^ , because of the different choices made both for 
the input EW parameters and parton distributions as 
detailed above. 



Much effort has been devoted to the study of Higgs 
production at the Tevatron at -^i = 2 TeV. These 
studies indicate that, given enough luminosity, a light 
Higgs boson can be discovered at the Tevatron using 
the associated production channels WH and ZH. In 
this report we present results of an analysis that 
incorporates as many of the backgrounds as possible 
at next-to-leading order for the WH channel. Whilst 
we use no detector simulation and do not attempt to 
include non-physics backgrounds, the results presented 
here can provide a normalization for more detailed 
studies. This is of importance since more detailed 
studies are often performed using shower Monte Carlo 
programs which can give misleading results for well 
separated jets. 

In particular, we will consider the light Higgs case 
(Af/z < 130 GeV) in the channel hbve^ . In 

addition to the usual cuts on rapidity and transverse 
momentum, 



\yh\Mh\ < 2, 

\ye\ < 2.5, 
\pI\,\pI\ > 15 GeV 
\pll\pt\ > 20 GeV 

we also impose isolation cuts, 

RbbJ R&b, Reb > ^■'^ I 



(29) 



(30) 



a s we ll as a cut on the scattering angle of the 66 system 
[ |l32|] ( the H iggs scattering angle) in the CoUins-Soper 
frame [ |l33tl , 



(31) 



I cos 6*^5 1 < 0.8 



Note that imposing the cut on cos Of^i requires knowl- 
edge of the longitudinal component of a neutrino 
momentum. Our results for the signal, backgrounds 
and significance are shown in Table ^, where we use 
e^g = 0.45 and integrate the cross-sections over a 
66 mass range appropriate for the Higgs mass under 
consideration. 



\Mh - A4b| < V2ctm, CTM - -IMh. 



(32) 



From this table, one can see that, even with a fairly 
restrictive set of cuts, the Wg* process in particular 
provides a challenging background. This is further 
emphasized in Figure ^, where the cross-sections for 
Mh = 110 GeV are presented in 5 GeV bins across 
the entire mf,^ spectrum. The signal, the two largest 
backgrounds and the sum of all the backgrounds 
including top quark production are plotted separately, 
as well as the totals with and without the Higgs signal. 
The sharp peak of the Higgs signal becomes only a 
small shoulder in the total distribution. 
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Table 3 

Signal, backgrounds (in fb) and significance for the VF-channel at y/s = 2 TeV 



i yj. v_Tt^ V 


Scalp 


100 


110 


120 


130 


bb) 


777, u 


8.8 


6.4 


4.2 


2.5 




(mjT. -A- 1 on Hf^V^ /9 


9^ 7 


99 7 




LO ~o 


W^Z{-^ bb) 


{mw + 100 GeV)/2 


6.7 


4.3 


2.0 


1.0 


bW+)t{^ feW^i^pJ 


100 GeV 


3.3 


3.7 


3.9 


3.9 




100 GeV 


0.3 


0.4 


0.5 


0.6 


W^*{t{-. bW+)b) 


100 GeV 


5.1 


5.8 


6.0 


6.0 


q't{-^ bW+) 


100 GeV 


0.3 


0.4 


0.5 


0.6 


Total B 




41.4 


37.3 


31.4 


27.6 


S/B 




0.21 


0.17 


0.13 


0.09 


s/Vb 




1.37 


1.05 


0.75 


0.48 



> 

<D 

CD 















PP ^ 


ve'bb with 


deteotor 


cuts 




and 


Mh = 110 GeV 






Wbb 




Total 














top 












25 50 


75 


100 


125 


150 
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Figure 9. Signal and backgrounds for WH. 'top' 
represents the sum of all the backgrounds including 
a top quark. 



5.4. Conclusions 

We have introduced the program MCFM, which 
calculates the rates for a number of standard model 
processes that are particularly relevant in Run II. 
These calculations are performed in fixed-order per- 
turbation theory, mainly at next-to-leading order in 
the strong coupling, and as such differ from other 
approaches such as parton shower Monte Carlos. As 
illustrations of the use of MCFM, we have presented 
total di-boson cross-sections and a primitive study 
(lacking detector effects and non-physics backgrounds) 
of WH production as a search for a light Higgs. Such 
calculations can be used to provide normalizations for 
more detailed studies in the future. 



6. Experimental handles on the backgrounds to 
new physics searches 

by Regina Demina 

6.1. Introduction 

Significant work has been done in the course of 
the SUSY/Higgs [ Q and Strong Dynamics [ ^ 
Workshops to understand the Tevatron discovery po- 
tential for new physics. Several promising signatures 
have been identified and the discovery reach has been 
estimated. In these studies, it was assumed that 
the systematic error on the signal and background 
normalization is similar in size to the statistical error, 
which is about 10%. Thus, the systematic error of each 
individual background process must be kept under 5%. 
Though it is probably a reasonable assumption, this 
will not happen automatically and dedicated studies 
are needed to achieve this goal. In this paper, we 
review the most important backgrounds to new physics 
and ways to estimate them in signal-depleted control 
samples. 

6.2. New physics signatures 

Associated vector boson and heavy flavor jets pro- 
duction is probably the most promising signature 
for new physics searches at the Tevatron. Standard 
Model (SM) Higgs boson J_ 6^ , Supersymmetry [ ^ , 
technicolor and topcolor [^| and even extra-dimension 
[ |6^ signatures may appear in these channels. 

Tables ^ and ||, show examples of new physics 
processes that can produce W + 2 iet and -I- 3 or 
more jet signatures. From the experimental point of 
view, a 'W is usually a high pT lepton accompanied 
by a significant missing energy (e.g. CDF Run I cuts 
are Prie, > 20 GeV/c, fJr > 20 GeV [0). In that 
sense, the supersymmetric partner of - xt " looks 
like a W, except its transverse mass will be inconsistent 
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# 


Process 


Model 


Special features 


1 


WH, H ^hb 


SM Higgs 


Resonance in Mf,^ 


2 




Technicolor 


Resonance in My^^i and 


3 


P?, ^ W^TT^,TT^ cb 


Technicolor 


Resonance in M^^^ and M^i 


4 


X1X2, Xi ^i^Xi, X2 bbxi 


SUSY 


Mriif^T) inconsistent with W 


5 


ti,t ^ bW,t ^ tx\,t ^ cxl 


SUSY 


Mriif^T) inconsistent with W 



Table 4 

Potential new physics signatures in the + 2 jet channel. From the experimental point of view, a 'W is a high 
PT lepton accompanied by significant missing energy (e.g., CDF cuts are pT(e,/i) > 20 GeV/c, fJr > 20 GeV). In 
that sense, Xi looks like a 'W." 



# 


Process 


Model 


Special features 


1 

2 
3 
4 
5 


tt,t-^bW,i-* 1x1,1^1x1 

itj^bxt 

it, i -> blv 
99, g^ti 
Z'{Vs,Tit)^ti 


SUSY 
SUSY 
SUSY 
SUSY 
Topcolor 


MTiifr) inconsistent with W 
Mt{^$t) inconsistent with W 
Mxiifr) inconsistent with W 
Mt(^-^t) inconsistent with W 
Resonance in Mtt 



Table 5 

Potential new physics signatures in the + 3 or more jet channel. 



with the W hypothesis, but this will become obvious 
only when significant statistics is accumulated. Some 
models predict special features, like resonance behavior 
in the bb invariant mass, while others do not. 

Table || presents new physics processes that can pro- 
duce Z + 2 jet signatures. Here we assume that the Z 
decays to a pair of leptons. Usually, a Z mass window 
cut is applied. In that sense, the supersymmetric 
partner of the Z ~ X2 ~ looks like a Z in only some 
regions of SUSY parameter space. If the Z decays to 
a pair of neutrinos, it produces missing energy. In this 
case, all the processes presented in Table || produce 
a fx + 2 jet signature. Table ^ shows additional 
new physics processes that result in a fx signature. 
As we see, these channels are very important for new 
physics searches, and the Standard Model backgrounds 
to these signatures must be thoroughly understood 
before any claims of discovery are made. 



Wbb Wbb 



1 . e or ^ with p-^> 20 GeV/c 

2. B-|,> 20 GcV 

3. 2 jets with E-j>15 GeV 

4. 2 b-tags 




Figure 10. Selection cuts and composition of the Wbb 
sample. 



6.3. Backgrounds to new physics 

The W{Z)bb signature was studied in the course 
of the SUSY/Higgs Workshop for the Higgs discovery 
potential estimate [ |6^. The ^T+heavy flavor (c/6) 
signature was studied in the CDF stop/sbottom search 
[ 1?^. We use these analyses as examples in our 
discussion. 



W{Z)bb signature 

Selection and sample composition 

The Wbb (Z(— > i'i>)bb) selection criteria and the re- 
sultant sample composition are summarized in Fig. ^ 
(Fig. |l^). The dominant contribution to both samples 
is QCD production of a vector boson accompanied by 
two 6-jets. 
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# 


Process 


Model 


Special features 


1 


ZH, H -fhb 


SM Higgs 


Resonance in Mj^i 


2 




Technicolor 


Resonance in M^cb and M^i 


3 




SUSY 


Mu inconsistent with Z, extra fJr 


4 


hh;b^bxlb^bxlxl^Uxl 


SUSY 


extra fx 



Table 6 

Potential new physics signatures in the Z + 2 jet channel. From the experimental point of view, a "Z" is two high 
Pt leptons, usually with a Z mass window cut. In that sense, X2 looks like a "Z" only in some regions of SUSY 
parameter space. 



# 


Process 


Model 


Special features 


1 


U,t^cx\ 




SUSY 


2 charm jets and Ipr 


1 


bb,b^b^. 




SUSY 


2 bottom jets and ]pT 


3 


LQ2LQ2,LQ2- 


cv 


Leptoquarks 


2 charm jets and IpT 


3 


LQ3LQ3,LQ3- 


bv 


Leptoquarks 


2 bottom jets and fx 



Table 7 

Potential new physics signatures in the IpT + 2 jet channel. This does not include processes complementary to those 
in Table ||, where a Z decays to a pair of neutrinos Z — » vv, thus producing missing energy. 



1 . E-j> 40 GeV 

2. 2jels withE-j>15 GcV 

3. min(AlKjet,E^))>0.5 

4. 2 b-tags 




ZH(llO) 



Figure 11. Selection cuts and composition of the Ipxbb 
sample. 



QCD W{Z)bb production. Experimental stud- 
ies of gluon splitting to heavy flavor. 

Diagrams of QCD associated production of W{Z) 
and heavy flavor jets are presented in Figure The 
leading contribution is Vl^(Z)+gluon production with 
subsequent gluon splitting to a 66 or cc pair, shown in 
Figure [l^(a). 

Though a next-to-leading-order calculation of the 
QCD Wbb production exists [ |7^, even the authors 
themselves recommend that it should be tested ex- 
perimentally. This is a particularly hard task in the 



presence of a potential signal contribution. In the case 
of the Higgs search, an invariant mass of two 6-jets 
could be used as an additional handle, since gluon 
splitting contributes mainly to the low part of the Mf,^ 
spectrum, while the Higgs is a resonance at 110-130 
GeV/c^. This is not the case for some other potential 
signal process, e.g. process 5 in Table ^. 

The probability for a gluon to split to two heavy 
flavor jets can be studied experimentally in different 
samples. The signal contamination becomes negligible, 
if the presence of a vector boson is not required. 

Three heavy flavor production mechanisms can be 
isolated - direct production, final state gluon splitting 
and initial state gluon splitting, also called flavor 
excitation. Diagrams of these processes are presented 
in Figure |l^. 

Though direct production is the lowest order pro- 
cess, it is responsible for the production of only ~ 20% 
of heavy flavor jets with energy above 20 GeV; about 
35% are produced by flavor excitation and 45% by 
gluon splitting. The relative contribution of different 
processes changes after 6-tagging is applied. Tagging is 
usually more efficient on directly produced jets, which 
tend to be back-to-back in the azimuthal plane. Heavy 
flavor quarks produced from gluon splitting are not 
well separated, and are often assigned to the same 
jet. Thus the relative contribution of gluon splitting 
to the double-tagged jet sample is quite low. Flavor 
excitation involves an initial state gluon splitting to 
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(a) 



(b) 



(c) 






Figure 12. Diagrams of QCD associated production of W{Z) and heavy flavor jets. 

Direct production Gluon splitting Flavor excitation 







Figure 13. Diagrams of QCD heavy flavor production. 



two heavy flavor quarks, one of which undergoes 
hard scattering. The other quark, being a part of 
the proton remnant, is often outside the detector 
acceptance. Thus, the contribution of flavor excitation 
to the double-tagged sample is significantly depleted. 
An analysis of the angular correlation between two 
heavy flavor tagged jets can be used to isolate the 
giuon-splitting component in heavy flavor production, 
as depicted in Fig. ^ 

Different methods can be used to tag heavy flavor 
jets. 

1. Impact parameter or secondary vertex tagging 
(JETPROB or SECVTX in CDF jargon) [ |7§ 
are the ones most commonly used. These 
samples have relatively high statistics. Using 
the same tagging method for the background 
and the signal sample eliminates the systematic 
uncertainty. The main disadvantage of these 
methods is the relatively low purity of these 
samples - contamination from c-jets and mistags 
is non-negligible. Usually, to get a stable fit, the 
h/c ratio has to be fixed to the one predicted 
by Monte Carlo, which is not without its own 
uncertainty. 

2. One of the heavy flavor jets is tagged by the 



presence of a high pT > 8 GeV/c lepton - 
electron or muon - and JETPROB or SECVTX 
tags another jet [ |7^. These samples have high 
statistics as well, but again suffer from charm 
and mistag contamination. Nonetheless, it is an 
interesting independent study. 

Both heavy ffavor jets can be tagged by leptons. 
In this case it is possible to go lower in lepton mo- 
mentum, usually > 3 GeV/c [Q. Compared 
to the first two cases, these studies probe a lower 
energy region, where the direct production mech- 
anism dominates. Thus not much information 
about gluon splitting probability can be gained. 

Study #2 can be modified to increase the purity 
by reconstructing exclusive or semi-exclusive fi- 
nal states in one of the jets: 

(a) Reconstructing a decay chain D* 

D^Ti, Ke(u)v can isolate the charm 

contribution [ . The presence of a high 
Pt lepton guarantees that the contribution 
from b ^ c decay is at the order of 10% 
or lower. Studying the angular correlation 
between D* jet and an impact parameter 
tagged jet isolates the gluon splitting to 
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Top, single top and diboson production 

Other backgrounds to new physics searches are top 
pair [ , single top [ Q and diboson [ ^ production. 
The theoretical predictions for these backgrounds are 
more reliable, because they do not involve gluon 
radiation and splitting, yet they still have to be tested 
experimentally. This is more or less a straightforward 
task for top pair and diboson production, where final 
states can be exclusively identified. It is less so for 
single top production, where the final state is exactly 
the same - Wbb - as in new physics channels in Table ^ 
Additional mass constraints, e.g. on the Wb mass can 
be used to isolate this process, but it is not at all 
obvious that adequate uncertainty can be reached for 
this channel. 

Missing energy and heavy flavor signatures 

Here, we summarize the selection criteria and com- 
position of the missing energy and heavy flavor sample 
used for top squark searches [ ^ : 



Figure 14. Distribution Ac/) between two 6-quarks from 
Monte Carlo events. 



charm contribution. This probability can 
then be applied to study #2 to extract the 
probability of gluon splitting to 6-quarks. 

(b) A similar trick can be used to isolate the b- 
contribution in the lepton sample. Promis- 
ing decay chains are [ |7^: B IvD* , 
D* D^TT, D° Kit or KSir; B -> lyD+, 
D+ Knir; B IvD", D" Ktt; and 
B Jf-^K. 

In Table |^, we present the number of events in each 
of the discussed samples collected in Run I and ex- 
pected in Run II and associated statistical uncertainty 
on the gluon splitting probability. 

Run I numbers are based on CDF results. In 
Run II, both CDF and D0 will have similar tracking 
and vertexing capabilities, thus these numbers are 
applicable to both detectors. Statistics in Run II is 
increased by a factor of 40, where 20 is gained from the 
luminosity increase and 2 from increased acceptance 
of the silicon microvertex detectors. The tagging 
efficiency increase is not taken into account. With 
these dedicated studies, the statistical uncertainty on 
the probability of gluon splitting to heavy flavor quarks 
can be significantly reduced in Run II, and will become 
adequate to the needs of new physics searches. 
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ww.wz.zz 



Figure 15. Selection cuts and composition of frc 
sample. 



More than 50% of the background is composed of 
W{-^ Tv) + 1 jet events. 

W + c production 

The leading order production process for W{—>- tv) + 
1 jet, where this jet is identified as charm, is sg — » Wc. 
The main uncertainty of the production rate for this 
process comes from the PDF of sea s quarks fs{x), 
which is measured by NuTeV [ ^ in the neutrino 
scattering process v^s ^ ^c. 

Figure |l6| (a) shows the distribution in x of sea 
s-quarks that contribute to sg — > Wc production at 
the Tevatron after the selection cuts from Figure fsl 
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# 


Sample 


N(Run I) 




-> QQ) Run I 


N(Run II) 


o{9 QQ) Run II 


1 
i 


Double tagged jets 






zU/o 


ZoUUU 


o.z/o 


2 


Muon+JETPROB tag 


2620 




16% 


104800 


2.5% 


3 




18000 




15% 


720000 


2.4% 


4 


B ^ D*£iy D"{^ Kit / KiTT)tuTT 


1700 




n.a. 


68000 


n.a. 


5 


B ^ D+{-^ KTnT)lv 


1900 




n.a. 


76000 


n.a. 


6 


B D°{-^ Kii)tv 


2700 




#4-#7 


108000 


#4-#7 


7 


B J/*i^(*) 


1300 




23% 


52000 


3.6% 



Table 8 

Data samples for heavy flavor production study. Numbers in samples ^1 and ^2 are double tags, while in samples 
=f^3-^7 numbers of exclusively reconstructed events are shown, without requiring a tag on the opposite side. 
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+ 1+ 

t u 11 1 
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b) 




Figure 16. (a) Distribution in x of sea s-quarks that 
contribute to sg — s- Wc production at the Tevatron, 
generated with PYTHIA 6.1+CTEQ4LO. Selection 
cuts have been applied, (b) Solid lines represent fs{x) 
and its uncertainty, as measured by NuTeV. It is 
compared to GRV94LO (dashed line) and CTEQ4L0 
(dot-dashed line) strange sea distributions. 



have been applied. Figure ^ (b) shows the fs{x) and 
its uncertainty measured by NuTeV. As we see, the 
region of NuTeV sensitivity is relevant for Tevatron 
studies. The overall uncertainty on fs{x) is 13.5%. 
fs (x) measured by NuTeV is in a good agreement with 
the results of CCFR [ |8l| , which has an uncertainty of 
10.5%. 

Since these uncertainties are dominated by the 
experimental statistics, we can hope that the combined 
result will have an uncertainty near 8.5%. The strange 
sea parton density function was also measured by the 



CHARM II [ H and CDHS [ || experiments. Combi- 
nation of results of all four experiments is certainly 
desirable, but non-trivial, since somewhat different 
techniques where used in each analysis. 

In Figure |l^(b), fs(x) measured by NuTeV is 
compared to the one provided by the GRV94LO [ 
H and CTEQ4L0 [ || PDF's, which are shown 
by dashed and dot-dashed lines, respectively. None 
of the PDF's provide an adequate description of the 
strange sea data. In the Tevatron search experiments, 
the systematic uncertainty on the background due 
to PDF's was typically estimated by switching from 
one PDF to another. In this case, the systematic 
uncertainty on the number of Wc events that pass our 
cuts is 36%. If instead of CTEQ4L0, fs{x) measured 
by NuTeV were used, the number of expected Wc 
events would go down by 30%. This is within the 
estimated uncertainty, but clearly the uncertainty has 
been overestimated. The correct uncertainty to use 
is 13.5%, or 8.5%, when the results of NuTeV and 
CCFR will be combined. This is a significant reduction 
compared to 36%, yet it is still not down to the 
desirable level of 5%. We can probably do better by 
studying Wc production when the W is identified by its 
decay to a muon or an electron. The expected number 
of events in the Wc, W £iy{£ = e, /i) channel is about 
2800, after applying the cuts listed in Figure |l5|, which 
corresponds to the statistical uncertainty for this back- 
ground of about 2%. The systematic uncertainty on 
missing energy and charm identification are common 
to the two channels, and the difference is in lepton vs. 
tau identification uncertainty, which can be expected 
to be below 5% with Run II statistics. 

QCD background 

The next dominant background in the ^tc channel is 
QCD multijet production, where missing energy comes 
from jet energy mis-measurement. This background 
is the toughest one to estimate, because it involves 
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multiple gluon radiation and splitting. Not only the 
overall rate, but also the angular correlation between 
jets may not be predicted reliably. To isolate this 
component, the usual trick is to apply all the signal 
selection cuts except tagging, subtract other known 
backgrounds and call the rest "QCD." The tagging 
probability derived from an independent jet sample 
is then applied to estimate the QCD contribution to 
the tagged sample [ |6|. One obvious drawback is 
that the heavy flavor fraction can change after the cuts 
are applied. Another is that the signal contribution is 
not always negligible even before tagging, and to some 
extent, it may be normalized away. 

Other backgrounds in the IfiTC channel. 

Other backgrounds in the channel are top pair 
production, di-bosons, Z+ jets and jets, where 
leptons were not identified. The discussion of these 
processes in Section 6.3 applies to the channel as 
well. 

7. Variable flavor number schemes for heavy 
quark electroproduction 

hy J. Smith 

Heavy quark production has been a major topic of 
investigation at hadron-hadron, electron-proton and 
electron-positron colliders. Here a review is given 
of some topics which are of interest primarily for 
electron-proton colliders. We concentrate on this 
reaction because a theoretical treatment can be based 
on the operator product expansion, and also because 
data are available for deep-inelastic charm production 
at HERA. How all this relates to Fermilab experiments 
will be discussed at the end. 

In QCD perturbation theory, one needs to introduce 
a renormalization scale and a mass factorization scale 
to perform calculations. We choose both equal to /i^, 
which will be a function of and the square of the 
mass of the charm quark m^. At small /i^, where 
kinematic effects due to quark masses are important, 
the best way to describe charm quark production is via 
heavy quark pair production from light quark u, d, s 
and gluon initial states. The mass m only appears in 
the heavy quark coefficient functions (or partonic cross 
sections) like Hff\z,Q'^,m'^,^j,'^), etc., [ ||]. Here 
the superscripts refer to their flavor decomposition and 
the order in perturbation theory, while the subscripts 
refer to the projection i = 2, L and the partonic 
initial state. The arguments refer to the partonic 
Bjorken variable z — Q'^/{s + Q^) and to the fact 
that these functions depend upon invariants and scales. 
The renormalization necessary to calculate these NLO 
expressions follows the CWZ method [ pS| . The symbol 



H refers to those coefficient functions which are derived 
from Feynman diagrams where the virtual photon 
couples to a heavy quark line. Analytic expressions for 
these functions are not known, but numerical fits are 
available in [ |89| . Asymptotic expressions in the limit 
^ are available in [ pG] . These contain terms 
like ln2(Q2/m2) and ln{Q^ /1^)Iti{Q^ / n^) multiplied 
by functions of z; they are manifestly singular in the 
limit that m — > 0. 

There are other heavy quark coefficient functions 



such as L. 



NS,(2) 
i,q 



z,Q^,TO^), which arise from tree di- 



agrams where the virtual photon attaches to the an 
initial state light quark line, so the heavy-quark is pair 
produced via virtual gluons. Analytic expressions for 
these functions are known for all z, and m? , which, 
in the limit ^ contain powers of ln(Q^/TO^) 
multiplied by functions of z. The three-flavor light 
mass MS parton densities can be defined in terms of 
matrix elements of operators and are now available 
in parton density sets. This is a fixed order per- 
turbation theory (FOPT) description of heavy quark 
production with three- flavor parton densities. Due 
to the work in [ the perturbation series is now 
known up to second order. In regions of moderate 
scales and invariants, this NLO description is well 
defined and can be combined with a fragmentation 
function to predict exclusive distributions [ |9^ for 
the outgoing charm meson, the anti-charm meson and 
the additional parton. This NLO massive charm 
approach agrees well with the recent ZJ-meson inclusive 
data in [ |9^ and [ The charm quark structure 

functions in this NLO description will be denoted 
if™(^,Q^m^n^ = 3). 

A different description, which should be more appro- 
priate for large scales where terms in are negligible, 
is to represent charm production by a parton density 
/c(a;,/i^), with a boundary condition that the density 
vanishes at small values of /.t^. Although at first sight 
these approaches appear to be completely different, 
they are in fact intimately related. It was shown in 
[ Q that the large terms in ln((5^/m^) which arise 
when ^ m^, can be resummed to all orders 
in perturbation theory. In this reference, all the 
two-loop corrections to the matrix elements of massive 
quark and massless gluon operators in the operator 
product expansion were calculated. These contain the 
same type of logarithms mentioned above multiplied 
by functions of z (which is the last Feynman inte- 
gration parameter). After operator renormalization 
and suitable reorganization of convolutions of the 
operator matrix elements (OME's) and the coefficient 
functions, the expressions for the infrared-safe charm 
quark structure functions Fi^c{x, Q^, m^. A) take on a 
simple form. After resummation, they are convolutions 
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of light-mass, four-flavor parton coefficient functions, 
commonly denoted by expressions like Cf'g'^\Q'^/fi'^) 
(available in [ [ with four-flavor light-parton 
densities, which also include a charm quark density 
fc(x,iJ?). Since the corrections to the OME's con- 
tain terms in hi{Q'^/m?) and ln(m^//x^) as well as 
non-logarithmic terms, it is simplest to work in the 
MS scheme with the scale jj?' = w? for < 
and ^2 = m2 Q^{1 - for > 

and discontinuous matching conditions on the flavor 
densities at /i^ = m^. Then all the logarithmic terms 
vanish at = ^"^ — m?' and the non-logarithmic 
terms in the OME's are absorbed into the boundary 
conditions on the charm density, the new four-flavor 
gluon density and the new light-flavor u,d,s densities. 
The latter are convolutions of the previous three-flavor 
densities with the OME's given in the Appendix of [ 

The above considerations lead to a precise descrip- 
tion through order of how, in the limit m — *■ 
0, to re-express the Ff^^'^^ {x,Q'^ .rri^) written in 
terms of convolutions of heavy quark coefficient func- 
tions with three-flavor light parton densities into a 
description in terms of four-flavor light-mass parton 
coeflacient functions convoluted with four-flavor par- 
ton densities. This procedure leads to the so-called 
zero-mass variable-flavor-number scheme (ZM-VFNS) 
for Fi^c{x,Q'^ 1 where the m dependent logarithms 
are absorbed into the new four-flavor densities. To 
implement this scheme, one has to be careful to use 
inclusive quantities which are coUinearly finite in the 
limit m ^ and A is an appropriate parameter 
which enables us to do this. In the expression for 
Fi,c there is a cancellation of terms in \T?{Q'^/m?') 
between the two-loop corrections to the light quark 
vertex function (the Sudakov form factor) and the 
convolution of the densities with the soft part of 
the L2,g-coefficient function. This is the reason for 
the split of Li^q into soft and hard parts, via the 
introduction of a constant A. Details and analytic 
results for LfO^^ LHARd ^^^.g available in [0. AU 
this analysis yielded and used the two-loop matching 
conditions on variable-flavor parton densities across 
flavor thresholds, which are special scales where one 
makes transitions from say a three-flavor massless 
parton scheme to a four-flavor massless parton scheme. 
The threshold is a choice of fi which has nothing 
to do with the actual kinematical heavy flavor pair 
production threshold at Q'^{x~^ - 1) = 4™^. In [ |4),[ 
H it was shown that the F^^^^^'^(x, Q^ m^, n/ = 3) 
tend numerically to the known asymptotic results in 
^asymP(^^q2^^2^„^ = 3), when » m^, which 

also equal the ZM-VFNS results. The last description 
is good for large (asymptotic) scales and contains 



a charm density fc{x,fi^) which satisfies a specific 
boundary condition at /z^ = m^. We denote the 
charm quark structure functions in this description by 

For moderate values of Q^, a third approach has 
been introduced to describe the charm components 
of Fi{x,Q'^). This is called a variable flavor number 
scheme (VFNS). A first discussion was given in [ p9|] , 
where a VFNS prescription called ACOT was given 
in lowest order only. A proof of factorization to all 
orders was recently given in [ IOC] for the total struc- 
ture functions Fi{x,Q'^), but the NLO expressions for 
Fi,c{x,Q^,m^) in this scheme were not provided. An 
NLO version of a VFNS scheme has been introduced 



97| and will be called the CSN scheme. A different 



approach, also generalized to all orders, was given 
in [ ||,[ |8), which is called the BMSN scheme. 
Finally another version of a VFNS was presented in 



101], which is called the TR scheme. The differences 



between the various schemes can be attributed to two 
ingredients entering the construction of a VFNS. The 
first one is the mass factorization procedure carried 
out before the large logarithms can be resummed. 
The second one is the matching condition imposed on 
the charm quark density, which has to vanish in the 
threshold region of the production process. All VFNS 
approaches require two sets of parton densities. One 
set contains three-flavor number densities whereas the 
second set contains four-flavor number densities. The 
sets have to satisfy the MS matching relations derived 
in J Appropriate four-flavor densities have been 
constructed in J ^ starting from the three-flavor LO 
and NLO sets of parton densities recently published in 
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Since the formulae for the heavy quark structure 
functions are available in J we only mention 

a few points here. The BMSN scheme avoids the 
introduction of any new coefficient functions other 
than those above. Since the asymptotic limits for 
3> of all the operator matrix elements and 
coefficient functions are known, we define (here Q 
refers to the heavy charm quark) 



F^^^^{x,Q\m',A,nj 



4) 



F, 



EXACT 



(x. 



)^m^A,n/ 



3) 



F, 



ASYMP 



(a;,Q2,m2, A,7i/ = 3) 



+ F[^^''{x,Q',m^,A,nf = A). (33) 

The scheme for F^q^ introduces a new heavy quark 
OME Aqq^^^(z, fi'^/m'^) J|l0| and coefficient functions 



Hf,Q^'^\z,Q^/m^) \ PI because it requires an incom- 



ing heavy quark Q, which did not appear in the NLO 
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Figure 17. The charm quark structure functions 
^EXACT(^^ = 3) (sohd hnc) F^f^ijif = 4), (dot- 
dashed hue) F^'^'^'^ijif ^ 4), (dashed hue) and 
Fl^^{nf = 4), (dotted hue) in NNLO for x = 0.005 
plotted as functions of Q^. 




Figure 18. The charm quark structure functions 
FE^acT(„^ = 3) (sohd hue) F^f^iuf = 4), (dot- 
dashed hue) F]^^^^{nf = 4), (dashed hue) and 
Pl^irtf = 4), (dotted Une) in NNLO for x = 0.005 
plotted as functions of Q^. 



corrections in [ |87| . The CSN coefficient functions are 
defined via the following equations. Up to second order 
we have 

■,CSN,SOFT,NS,(2) / » \ _ ' 



NS,(1) ( _^VIRT,NS,(2)|. Q^sj 



2 n2s 



-C™'^^^(a4-^), (34) 

with the virtual term the second order Sudakov form 
factor. The other CSN coefficient functions are defined 
by equations like (we only give one of the longitudinal 
terms for illustration) 

^CSN,S,(l)/Q^ ^ 
2 



r. n \ T ^an \ o ) r o 1 t 



(35) 



with C 



CSN,NS,(0) 



4toVQ2^ The CSN and BMSN 



schemes are designed to have the following two prop- 
erties. First of all, suppressing unimportant labels, 

^;T("/ = 4) = i^TX-/ = 4) 

= ^'q^''^(^/ = 3) for Q2<^2 (3g) 

Since /q(to^)^^^° 7^ (see [||l) this condition can be 
only satisfied when we truncate the perturbation series 



at the same order. The second requirement is that 



hm i^P,^SN 



(n/ = 4)- hm ^;^^(n/=4) 



= hm i^^Mn/ = 4). (37) 

The only differences between the two schemes arises 
from terms in m? so they may not be equal just above 
— m? . This turns out to be the case for the 
longitudinal structure function, which is more sensitive 
to mass effects. 

Figure |l^ shows NNLO results for the (jg_ 
pendence of Ff^^^^{nf = 3), F^^'^{nf = 4), 
F2^^^{nf = 4)', and F^^^{nf = 4)'at x = 0.005. 
Note that the results satisfy the requirements in 
Eqs. (36) and (37). The ZM-VFNS description is 
poor at small . Figure |l8| shows the results for 
^L™("/ = 3), F^fiuf = 4), Ff^^^iuf 4), 
and Fl^^{nf = 4) at x = 0.005. We see that the 
CSN result is negative and therefore unphysical for 
2.5 < < 6 (GeV/c)2 which is due to the term in 
Aw? IQ^ and the subtraction in Eq. (35). 

One way this research work is of relevance to Fermi- 
lab experiments is that it produces more precise ZM- 
VFNS parton densities. Such densities are used ex- 
tensively to predict cross sections at high energies, for 
example for single top quarks. Therefore the previous 
work on four-flavor parton densities has been extended 
in [ 105 1 to incorporate the two- loop discontinuous 
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matching conditions across the bottom flavor threshold 
at fj, = rrib and provided a set of five-flavor densities, 
which contains a bottom quark density fi,{x,fjl^). The 
differences between the five-flavor densities and those 
and [ p5[ are also discussed. Results for 



106] 



deep-inelastic electroproduction of bottom quarks will 



be presented in [ 107 



8. The Underlying Event in Hard Scattering 
Processes 

by Rick Field and David Stuart 

8.1. Introduction 

The total proton-antiproton cross section is the sum 
of the elastic and inelastic cross sections. The inelastic 
cross section consists of a single-diffractive, double- 
diffractive, and a "hard core" component, where the 
"hard core" is everything else. "Hard core" does not 
necessarily imply "hard scattering." A "hard scatter- 
ing" collision, such as that illustrated in Fig. |l^(a), is 
one in which a "hard" (i.e. large transverse momen- 
tum) 2-to-2 parton-parton subprocess has occurred. 
"Soft" hard core collisions correspond to events in 
which no "hard" interaction has occurred. When there 
is no large pt subprocess in the collision, one is not 
probing short distances and it probably does not make 
any sense to talk about partons. The QCD "hard 
scattering" cross section grows with increasing collider 
energy and becomes a larger and larger fraction of the 
total inelastic cross section. In this analysis, we used 
the CDF Min-Bias trigger data sample in conjunction 
with the JET20 trigger data sample to study the 
growth and development of "charged particle jets" 
from pxijet) — 0.5 to 50GeV. We compared several 
"local" jet observables with the QCD "hard scattering" 
Monte-Carlo models of HERWIG[ §, IS A JET [ ||, 
and PYTHIA[|. 

A "hard scattering" event, like that illustrated in 
Fig. ^(a) consists of large-pr outgoing hadrons that 
originate from the large-pr partons {i.e., outgoing hard 
scattering "jets") and also hadrons that originate from 
the break-up of the proton and antiproton (i.e., the 
"beam-beam remnants"). The "underlying event" is 
an interesting object that is not very well understood. 
In addition to beam-beam remnants, it may contain 
hadrons resulting from initial-state radiation. Also, it 
is possible that multiple parton scattering occurs in 
hadron-hadron collisions as illustrated in Fig. |l9|(b). 
This is a controversial issue, but the underlying event 
might also contains hadrons that originate from mul- 
tiple parton interactions. PYTHIA, for example, uses 
multiple parton interactions as a way to enhance the 
activity of the underlying event [ |^ . 

In this analysis, we studied a variety of "global" 



observables to probe the growth and structure of 
the underlying event. We find that the underlying 
"hard scattering" event is not the same as a "soft" pp 
collision. For the same available energy, the underlying 
event in a hard scattering is considerably more active 
{i.e., higher charged particle density and more px) 
than a "soft" collision. This is not surprising since 
a violent hard scattering has occurred! We find that 
none of the QCD Monte-Carlo models {with their 
default parameters) describe correctly all the properties 
of the underlying event. 

8.2. Data Selection and Monte-Carlo Models 
(1) Data Selection 

The CDF detector, described in detail in Ref. [ 
122], measures the trajectories and transverse mo- 



menta, Pt, of charged particles in the pseudorapidity 
region \ri\ < 1.1 with the central tracking chamber 
(CTC), silicon vertex detector (SVX), and vertex time 
projection chamber (VTX), which are immersed in a 
1.4 T solenoidal magnetic field. In this analysis we 
consider only charged particles measured in the central 
tracking chamber (CTC) and use the two trigger 
sets of data listed in Table The minimum bias 
(min-bias) data were selected by requiring that at least 
one particle interacted with the forward beam-beam 
counter BBC (3.4 < 77 < 5.9) and/or the backward 
BBC (-5.9 < rj < -3.4). The min-bias trigger 
selects predominately the "hard core" component of 
the inelastic cross section. 

Charged particle tracks are found with high effi- 
ciency as long as the density of particles is not high. 
To remain in a region of high efficiency, we consider 
only charged particles with pT > 0.5 GeV and \ri\ < 1. 
The observed tracks include some fake tracks that 
result from secondary interactions between primary 
particles, including neutral particles, and the detector 
material. There are also particles originating from 
other pp collisions. To reduce the contribution from 
these sources, we consider only tracks which point 
to the primary interaction vertex within 2 cm along 
the beam direction and 1 cm transverse to the beam 
direction. Detector simulations indicate that this 
impact parameter cut is very efficient and that the 
number of fake tracks is about 3.5% when a 1 cm 
impact parameter cut is applied in conjunction with 
a 2 cm vertex cut. Without the impact parameter cut 
the number of fake tracks is approximately 9%. 

This dependence of the number of fake tracks on 
the CTC impact parameter cut provides a method 
of estimating systematic uncertainties due to fakes. 
Every data point P on every plot in this analysis was 
determined three times by using a 2 cm vertex cut 
in conjunction with three different CTC do cuts; a 
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Figure 19. (a) Illustration of a pp collision in which a "hard" 2-to-2 parton scattering with transverse momentum, 
PT(hard), has occurred. The resulting event contains particles that originate from the two outgoing partons (plus 
final-state radiation) and particles that come from the breakup of the p and p {i.e., "beam-beam remnants"). 
The "underlying event" consists of the beam-beam remnants plus initial-state radiation; (b) Illustration of a pp 
collision in which a multiple parton interaction has occurred. In addition to the "hard" 2-to-2 parton scattering 
with transverse momentum, p7^(hard), there is an additional "semi- hard" parton-parton scattering that contributes 
particles to the "underlying event." For PYTHIA, we include the contributions from multiple parton scattering in 
the beam-beam remnant component. 



Table 9 

Data sets and selection criterion used in this analysis. 



CDF Data Set 


Trigger 


Events 


Selection 


Min-Bias 


Min-Bias Trigger 


626,966 


zero or one vertex in \z\ < 100 cm 
\zc - Zy\ <2 cm, CTC do\ < 1 cm 
p^r^'"'^ > 0.5 GeV, |?7| < 1 


JET20 


Calorimeter Tower cluster 
with Et > 20 GeV 


78,682 


zero or one vertex in \z\ < 100 cm 
\zc - Zt, < 2 cm, CTC do\ < 1 cm 
p^-^'^'^ > 0.5 GeV, \r]\ < 1 
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1 cm CTC do cut (P), a 0.5 cm CTC da cut (Pi), 
and no CTC do cut (P2)- The 1 cm cut determined 
the value of the data point, P, and the difference 
between the 0.5 cm cut value and no cut value of 
the data point determined the systematic error of the 
data point as follows: sys-error = PjPj — Pi|/Pi This 
systematic error was then added in quadrature with 
the statistical error. We do not correct the data for the 
CTC track-finding efficiency. Instead, the theoretical 
Monte-Carlo model predictions are corrected. 

(2) QCD "Hard Scattering" Monte-Carlo 
Models 

The "hard" scattering QCD Monte-Carlo models 
used in this analysis are listed in Table The 
QCD perturbative 2-to-2 parton-parton differential 
cross section diverges as the pr of the scattering, p^*''^, 
goes to zero (see Fig. |l^). One must set a minimum 
p^'^'^ large enough so that the resulting cross section is 
not larger that the total "hard core" inelastic cross 
section, and also large enough to ensure that QCD 
perturbation theory is applicable. In this analysis we 
take p^^'^'i > 3 GeV. 

Each of the QCD Monte-Carlo models handle the 
"beam-beam remnants" in a similar fashion. A hard 
scattering event is basically the superposition of a hard 
parton-parton interaction on top of a "soft" collision. 
HERWIG[ I assumes that the underlying event is 
a soft collision between the two "beam clusters." 
ISAJET[||] uses a model similar to the one it uses for 
soft "min-bias" events {i.e., "cut Pomeron"), but with 
different parameters, to describe the underlying beam- 
beam remnants. PYTHIA[ ^ assumes that each in- 
coming beam hadron leaves behind a "beam remnant," 
which do not radiate initial state radiation, and simply 
sail through unaffected by the hard process. However, 
unlike HERWIG and ISAJET, PYTHIA also uses 
multiple parton interactions to enhance the activity of 



the underlying event as illustrated in Fig. 19 

In this analysis we examine two versions of 
PYTHIA, PYTHIA 6.115 and PYTHIA 6.125 both 
with the default values for all the parameters. The 
default values of the parameters are different in version 
6.115 and 6.125. In particular, the effective mini- 
mum pt for multiple parton interactions, PARP(81), 
changed from 1.4 GeV in version 6.115 to 1.9 GeV 
in version 6.125. Increasing this cut-off decreases 
the multiple parton interaction cross section which 
reduces the amount of multiple parton scattering. 
For completeness, we also consider PYTHIA with no 
muhiple parton scattering (MSTP(81)=0). 

Since ISAJET employs "independent fragmenta- 
tion" it is possible to trace particles back to their origin 
and divide them into three categories: particles that 



arise from the break-up of the beam and target (beam- 
beam remnants), particles that arise from initial-state 
radiation, and particles that result from the outgoing 
hard scattering jets plus final-state radiation. The 
"hard scattering component" consists of the particles 
that arise from the outgoing hard scattering jets plus 
initial and final-state radiation (sum of the last two 
categories). Particles from the first two categories 
(beam-beam remnants plus initial-state radiation) are 
normally what is referred to as the underlying event 
(see Fig. |l^). Of course, these categories are not 
directly observable experimentally. Nevertheless, it 
is instructive to examine how particles from various 
origins affect the experimental observables. 

Since HERWIG and PYTHIA do not use indepen- 
dent fragmentation, it is not possible to distinguish 
particles that arise from initial-state radiation from 
those that arise from final-state radiation, but we can 
identify the beam-beam remnants. When, for example, 
a color string breaks into hadrons it is not possible 
to say which of the two partons producing the string 
was the parent. For HERWIG and PYTHIA, we 
divide particles into two categories: particles that arise 
from the break-up of the beam and target (beam-beam 
remnants), and particles that result from the outgoing 
hard scattering jets plus initial and final-state radi- 
ation (hard scattering component). For PYTHIA, 
we include particles that arise from multiple parton 
interactions in the beam-beam remnant component. 

(3) Method of Comparing Theory with Data 

Our philosophy in comparing the theory with data in 
this analysis is to select a region where the data is very 
"clean." The CTC efficiency can vary substantially 
for very low px tracks and in dense high pT jets. 
To avoid this we have considered only the region 
Pt > 0.5 GeV and \ri\ < 1 where the CTC efficiency 
is high and stable (estimated to be 92% efficient) 
and we restrict ourselves to jets less than 50 GeV. 
The data presented here are uncorrected. Instead 
the theoretical Monte-Carlo predictions are corrected 
for the track finding efficiency and have an error 
(statistical plus systematic) of about 5%. The errors 
on the (uncorrected) data include both statistical and 
correlated systematic uncertainties. 

In comparing the QCD "hard scattering" Monte- 
Carlo models with the data, we require that the Monte- 
Carlo events satisfy the CDF min-bias trigger and we 
apply an 8% correction for the CTC track finding 
efficiency. The corrections are small. On the average, 
8 out of every 100 charged particles predicted by the 
theory are removed from consideration. Requiring 
the theory to satisfy the min-bias trigger is important 
when comparing with the Min- Bias data, but does not 
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Table 10 

Theoretical QCD "hard" scattering Monte-Carlo models studied in this analysis. In all cases we take pT(hard) > 
3GeV. 



Monte-Carlo Model 


Subprocess 


Comments 


HERWIG 5.9 


QCD 2-to-2 parton scattering 
IPROC = 1500 


Default values for all parameters 


ISAJET 7.32 


QCD 2-to-2 parton scattering 
TWOJET 


Default values for all parameters 


PYTHIA 6.115 


QCD 2-to-2 parton scattering 
MSEL = 1 


Default values for all parameters: 
PARP(81) = 1.4 


PYTHIA 6.125 


QCD 2-to-2 parton scattering 
MSEL = 1 


Default values for all parameters: 
PARP(81) = 1.9 


PYTHIA No MS 


QCD 2-to-2 parton scattering 
MSEL = 1 


Default values for all parameters: 

MSTP(81) = 



matter when comparing with the JET20 data since 
essentially all high pr jet events satisfy the min-bias 
trigger. 

8.3. The Evolution of Charge Particle "Jets" 
from 0.5 to 50 GeV 

We define charged particle "jets" and examine the 
evolution of these "jets" from p!^^ = 0.5 to 50 GeV. As 
illustrated in Fig. "jets" are defined as "circular 
regions" {R = 0.7) in ij-cj) space and contain charged 
particles from the underlying event as well as particles 
which originate from the fragmentation of high pT 
outgoing partons (see Fig. |l^). Also, every charged 
particle in the event is assigned to a "jet," with the 
possibility that some "jets" might consist of just one 
charged particle. We adapt a very simple jet definition 
since we will be dealing with "jets" that consist of 
only a few low px charged particles. The standard 
jet algorithm based on calorimeter clustering is not 
applicable at low px ■ 

(1) Jet Definition (charged particles) 

We define jets as circular regions in rj-cf) space with 
"distance" R = ^(ATy)^ + (A</«)2. Our jet algorithm 
is as follows: 

• Order all charged particles according to their pT- 

• Start with the highest pt particle and include in 
the "jet" all particles within R = 0.7. 

• Go to the next highest px particle [not already 
included in a "jet") and add to the "jet" all 
particles (not already included in a "jet") within 
R = 0.7. 

• Continue until all particles are in a "jet." 




Figure 20. Illustration of an event with six charged particles 
{Pt > 0.5 GeV and \r]\ < 1) and five charged "jets" (circular 
regions in ri-ip space with i? = 0.7). 



We consider all charged particles {pr > 0.5 GeV and 
I77I < 1) and allow the jet radius to extend outside 
1 77 1 < 1. Fig. ^ illustrates an event with six charged 
particles and five jets. We define the pT of the jet to 
be the scalar px sum of all the particles within the jet 
(i.e., it is simply the scalar sum within the circular 
region) . 

We know that the simple charged particle jet def- 
inition used here is not theoretically favored since if 
applied at the parton level it is not infrared safe. Of 
course, all jet definitions (and in fact all observables) 
are infrared safe at the hadron level. We have done 
a detailed study comparing the nave jet definition 
used here with a variety of more sophisticated charge 
particle jet definitions. This analysis will be presented 
in a future publication. Some of the observables 
presented here do, of course, depend on one's definition 
of a jet and it is important to apply the same definition 
to both the theory and data. 
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(2) Charged Jet Multiplicity versus pT(jet#l) 

Fig. ^ shows the average number of charged parti- 
cles {pt > 0.5 GeV and \r]\ < 1) within jet#l {leading 
charged jet) as a function of pT(jet#l). The solid 
points are Min-Bias data and the open points are the 
JET20 data. The JET20 data connect smoothly to 
the Min-Bias data and allow us to study observables 
over the range 0.5 < pT(jet#l) < 50 GeV. There 
is a small overlap region where the Min-Bias and 
JET20 data coincide. The errors on the data include 
both statistical and correlated systematic uncertain- 
ties, however, the data have not been corrected for 
efficiency. Fig. ^ shows a sharp rise in the leading 
charged jet multiplicity at low pT(jet#l) and then a 
flattening out and a gradual rise at high pT(jet#l). 
The data are compared with the QCD "hard scat- 
tering" Monte-Carlo predictions of HERWIG 5.9, 
ISAJET 7.32, and PYTHIA 6.115. The theory curves 
are corrected for the track finding efficiency and have 
an error {statistical plus systematic) of around 5%. 



defined in two ways: size according to particle number 
and size according to px- The first corresponds to the 
radius in rj-cf) space that contains 80% of the charged 
particles in the jet, and the second corresponds to 
the radius in -q-cj) space that contains 80% of the jet 
Pt- The data on the average "jet size" of the leading 
charge particle jet are compared with the QCD "hard 
scattering" Monte-Carlo predictions of HERWIG 5.9, 
ISAJET 7.32, and PYTHIA 6.115 in Fig. ||. A 
leading 20 GeV charged jet has 80% of its charged 
particles contained, on the average, within a radius in 
r]-(j) space of about 0.33, and 80% of its pr contained, 
on the average, within a radius of about 0.20. Fig. H] 
clearly illustrates the "hot core" of jets. The radius 
containing 80% of the pr is smaller than the radius 
that contains 80% of the particles. Furthermore, the 
radius containing 80% of the pr decreases as the overall 
Pt of the jet increases due to limited momentum 
perpendicular to the jet direction. 
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Figure 21. The average number of charged particles 
{pt > 0.5 GeV and \r]\ < 1) within the leading charged 
jet {R = 0.7) as a function of the pt of the leading 
charged jet. The solid (open) points are Min-Bias 
(JET20) data. The errors on the {uncorrected) data 
include both statistical and correlated systematic un- 
certainties. The QCD "hard scattering" theory curves 
(HERWIG 5.9, ISAJET 7.32, PYTHIA 6.115) are 
corrected for the track finding efficiency and have an 
error {statistical plus systematic) of around 5%. 



(3) Charged Jet "Size" versus pT(jet#l) 

Although the charged particle jets are defined as 
circular regions in rj-cj) space with R = 0.7, this is 
not the "size" of the jet. The "size" of a jet can be 



Figure 22. The average radius in rj-cf) space containing 
80% of the charged particles (and 80% of the charged 
Pt) as a function of the pT of the leading charged jet. 
The errors on the {uncorrected) data include both sta- 
tistical and correlated systematic uncertainties. The 
QCD "hard scattering" theory curves (HERWIG 5.9, 
ISAJET 7.32, PYTHIA 6.115) are corrected for the 
track finding efficiency and have an error {statistical 
plus systematic) of around 5%. 



8.4. The Overall Event Structure as a Function 

ofpT(jet#l) 

In the previous section, we studied "local" leading 
jets observables. The QCD Monte-Carlo models did 
not have to describe correctly the entire event in order 
to fit the observable. They only had to describe 
correctly the properties of the leading charge particle 
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jet, and all the models fit the data fairly well {although 
not perfectly) . Nowwe will study "global" observables, 
where to fit the observable the QCD Monte-Carlo 
models will have to describe correctly the entire event 
structure. 

(1) Overall Charged Multiplicity versus 

PT(jet#l) 

Figure ^ shows the average number of charged 
particles in the event with pT > 0.5 GeV and |?7| < 1 
(including jet#l) as a function of pT(jet#l) (leading 
charged jet) for the Min-Bias and JET20 data. Again 
the JET20 data connect smoothly to the Min-Bias 
data, and there is a small overlap region where the 
Min-Bias and JET20 data coincide. Figure ^ shows 
a sharp rise in the overall charged multiplicity at low 
PT(jet#l) and then a flattening out and a gradual rise 
at high pT(jet#l) similar to Fig. We would like to 
investigate where these charged particles are located 
relative to the direction of the leading charged particle 
jet. 
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Figure 23. The average total number of charged 
particles in the event (pr > 0.5 GeV and \ri\ < 1 
including jct#l) as a function of the pr of the leading 
charged jet. The solid (open) points are the Min- Bias 
(JET20) data. The errors on the (uncorrected) data 
include both statistical and correlated systematic un- 
certainties. The QCD "hard scattering" theory curves 
(HERWIG 5.9, ISAJET 7.32, PYTHIA 6.115) are 
corrected for the track finding efficiency and have an 
error (statistical plus systematic) of around 5%. 



(2) Correlations in relative to pr(jet#l) 

As illustrated in Fig. ^ the angle A(j> is defined 
to be the relative azimuthal angle between charged 



particles and the direction of the leading charged 
particle jet. We label the region |0 — 0jct#i| < 60° 
as "toward" jet#l and the region |0— 0jct#i| > 120° 
is as "away" from jet#l. The "transverse" to jet#l 
region is defined by 60° < |0 - 0jot#i| < 120°. Each 
region, "toward," "transverse," and "away" covers the 
same range |A?7| x \A(t)\ = 2 x 120°. The "toward" 
region includes the particles from jet#l as well as a 
few particles from the underlying event. As we will 
see, the "transverse" region is very sensitive to the 
underlying event. The "away" region is a mixture 
of the underlying event and the "away-side" hard 
scattering jet. 



Jet #1 Direction 




Figure 24. Illustration of correlations in azimuthal 
angle A(j> relative to the direction of the leading 
charged jet in the event, jet#l. The angle A(j> = 
(p — 0jct#i is the relative azimuthal angle between 
charged particles and the direction of jet#l. The 
region \A(f)\ < 60° is referred to as "toward" jet#l 
(includes particles in jet#l) and the region \A(j)\ > 
120° is called "away" from jet#l. The "transverse" to 
jet#l region is defined by 60° < \A(j)\ < 120°. Each 
region, "toward," "transverse," and "away" covers the 
same range |A?7| x \A<j)\ = 2 x 120°. 



Figure |25| shows the data on the average number 
of charged particles (pr > 0.5 GeV and \r]\ < 1) as a 
function of pT(jet#l) for the three regions. Each point 
corresponds to the "toward," "transverse," or "away" 
(Nchg) in a 1 GeV bin. The solid points are Min-Bias 
data and the open points are JET20 data. The data 
in Fig. |2^ define the average event "shape." For 
example, for an "average" pp collider event at 1.8 TeV 
with pT(jet#l) = 20 GeV there are, on the average, 
8.7 charged particles "toward" jet#l (including the 
particles in jet#l), 2.5 "transverse" to jet#l, and 4.9 
"away" fromjet#l. 
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Figure 25. The average number of "toward" (|A(/)| < 
60°), "transverse" (60° < |A0| < 120°), and "away" 
{\A(j)\ > 120°) charged particles {pr > 0.5 GeV and 
|?7| < 1 including jet#l) as a function of the pt 
of the leading charged jet. Each point corresponds 
to the {Nchg) in a 1 GeV bin. The sohd (open) 
points are the Min-Bias (JET20) data. The errors 
on the (uncorrected) data include both statistical and 
correlated systematic uncertainties. The "toward," 
"transverse," and "away" regions are defined in Fig. 



Figure ^ shows the data on the average scalar pT 
sum of charged particles {px > 0.5 GeV and \r]\ < 1) 
as a function of pT(jct#l) for the three regions. Here 
each point corresponds to the "toward," "transverse," 
or "away" (pTsum) in a 1 GeV bin. In Fig. |2^, data on 
(Nchg) as a function of p7-(jet7^1) for the three regions 
are compared with the QCD "hard scattering" Monte- 
Carlo predictions of HERWIG 5.9, ISA JET 7.32, and 
PYTHIA 6.115. The QCD Monte-Carlo models agree 
qualitatively (but not precisely) with the data. We will 
now examine more closely these three regions. 

(3) The "Toward" and "Away" Region versus 

PT(jet#l) 

Figure |2^ shows the data from Fig. ^ on the 
average number of "toward" region charged particles 
compared with the QCD "hard scattering" Monte- 
Carlo predictions of HERWIG 5.9, ISAJET 7.32, 
and PYTHIA 6.115. This plot is very similar to 
the average number of charged particles within the 
leading jet shown in Fig. |p. At pT(jet#l) = 20 GeV, 
the "toward" region contains, on the average, about 
8.7 charged particles with about 6.9 of these charged 
particles belonging to jet#l. As expected, the toward 
region is dominated by the leading jet. This is seen 
clearly in Fig. ^ where the predictions of ISAJET for 
the "toward" region are divided into three categories: 
charged particles that arise from the break-up of the 



Figure 26. The average scalar pr sum of "toward" 
(|A</)| < 60°), "transverse" (60° < |A</)| < 120°), and 
"away" {\A(t)\ > 120°) charged particles {pr > 0.5 GeV 
and \ri\ < 1 including jet#l) as a function of the pr 
of the leading charged jet. Each point corresponds 
to the (PTsum) in a 1 GeV bin. The solid (open) 
points are the Min-Bias (JET20) data. The errors 
on the (uncorrected) data include both statistical and 
correlated systematic uncertainties. The "toward," 
"transverse," and "away" regions are defined in Fig. 
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Figure 27. Data from Fig. |25| on the average number 
of "toward" (\A(f>\ < 60°), "transverse" (60° < |A0| < 
120°), and "away" (|A0| > 120°) charged particles 
(pt > 0.5 GeV and |?7| < 1 including jet#l) as a 
function of the pr of the leading charged jet compared 
to QCD "hard scattering" Monte-Carlo predictions of 
HERWIG 5.9, ISAJET 7.32, and PYTHIA 6.115. 
The errors on the (uncorrected) data include both sta- 
tistical and correlated systematic uncertainties. The 
theory curves are corrected for the track finding effi- 
ciency and have an error (statistical plus systematic) of 
around 5%. 
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beam and target {beam-beam remnants), charged parti- 
cles that arise from initial-state radiation, and charged 
particles that result from the outgoing jets plus final- 
state radiation. For pT(jet#l) values below 5GeV 
the "toward" region charged multiplicity arises mostly 
from the beam-beam remnants, but as pT(jet#l) 
increases the contribution from the outgoing jets plus 
final-state radiation quickly begins to dominate. The 
bump in the beam-beam remnant contribution at low 
PT(jet#l) is caused by leading jets composed almost 
entirely from the remnants. 
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Figure 28. Data from Fig. ^ on the average number 
of charged particles {pr > 0.5 GeV and |?7| < 1) as 
a function of pT(jet#l) {leading charged jet) for the 
"toward" region " defined in Fig. 24 compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of HERWIG 5.9, ISAJET 7.32, and PYTHIA 6.115. 
Each point corresponds to the "toward" {Nchg) in 
a 1 GeV bin. The errors on the {uncorrected) data 
include both statistical and correlated systematic un- 
certainties. The theory curves are corrected for the 
track finding efficiency and have an error {statistical 
plus systematic) of around 5%. 



Fig. |30| shows the data from Fig. g5| on the aver- 
age number of "away" region charged particles com- 
pared with the QCD "hard scattering" Monte-Carlo 
predictions of HERWIG 5.9, ISAJET 7.32, and 
PYTHIA 6.115. In Fig. the data from Fig. || 
on the average scalar px sum in the "away" region is 
compared to the QCD Monte-Carlo predictions. The 
"away" region is a mixture of the underlying event 
and the "away-side" outgoing "hard scattering" jet. 
This can be seen in Fig. |3^ where the predictions 
of ISAJET for the "away" region are divided into 
three categories: beam-beam remnants, initial-state 
radiation, and outgoing jets plus final-state radiation. 



Figure 29. Data from Fig. ^ on the average number 
of charged particles {pt > 0.5 GeV and I77I < 1) as 
a function of pT(jet#l) {leading charged jet) for the 
"toward" region defined in Fig. |2j compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of ISAJET 7.32. The predictions of ISAJETare 
divided into three categories: charged particles that 
arise from the break-up of the beam and target 
{beam-beam remnants), charged particles that arise 
from initial-state radiation, and charged particles that 
result from the outgoing jets plus final-state radiation 
(see Fig. |l^). The errors on the {uncorrected) data 
include both statistical and correlated systematic un- 
certainties. The theory curves are corrected for the 
track finding efficiency and have an error {statistical 
plus systematic) of around 5%. 



Here the underlying event plays a more important role 
since the "away-side" outgoing "hard scattering" jet 
is sometimes outside the regions |?7| < 1. For the 
"toward" region the contribution from the outgoing 
jets plus final state-radiation dominates for pr(jet#l) 
values above about 5 GeV, whereas for the "away" 
region this does not occur until around 20 GeV. 

Both the "toward" and "away" regions are described 
fairly well by the QCD "hard scattering" Monte-Carlo 
models. These regions are dominated by the outgoing 
"hard scattering" jets and as we saw in Section C the 
Monte-Carlo models describe the leading outgoing jets 
fairly accurately. We will now study the "transverse" 
region which is dominated by the underlying event. 

8.5. The "Transverse" Region and the Under- 
lying Event 

Fig. shows that there is a lot of activity in the 
"transverse" region. If we suppose that the "trans- 
verse" multiplicity is uniform in azimuthal angle </> and 
pseudo-rapidity 77, the observed 2.3 charged particles at 
PT(jet#l) = 20 GeV translates to 3.8 charged particles 
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Figure 30. Data from Fig. ^on the average number of 
charged particles {pr > 0.5 GeV and |?7| < 1) as a func- 
tion of pT(jet#l) {leading charged jet) for the "away" 
region defined in Fig. ^ compared with the QCD 
"hard scattering" Monte-Carlo predictions of HER- 
WIG 5.9, ISAJET 7.32, and PYTHIA 6.115. The 
errors on the (uncorrected) data include both statistical 
and correlated systematic uncertainties. The theory 
curves are corrected for the track finding efficiency and 
have an error (statistical plus systematic) of around 5%. 
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Figure 31. Data from Fig. on the average scalar pr 
sum of charged particles {{pr > 0.5 GeV and \r]\ < 1) 
as a function of pT(jet#l) (leading charged jet) for 
the "away" region defined in Fig. |2 compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of HERWIG 5.9, ISAJET 7.32, and PYTHIA 6.115. 
The errors on the (uncorrected) data include both sta- 
tistical and correlated systematic uncertainties. The 
theory curves are corrected for the track finding effi- 
ciency and have an error (statistical plus systematic) of 
around 5%. 



Figure 32. Data from Fig. ^ on the average number 
of charged particles (px > 0.5 GeV and \r]\ < 1) 
as a function of prOet^l) (leading charged jet) for 
the "away" region defined in Fig. |2j compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of ISAJET 7.32. The predictions of ISAJETare 
divided into three categories: charged particles that 
arise from the break-up of the beam and target 
(beam-beam remnants), charged particles that arise 
from initial-state radiation, and charged particles that 
result from the outgoing jets plus final-state radiation 
(see Fig. |l^). The errors on the (uncorrected) data 
include both statistical and correlated systematic un- 
certainties. The theory curves are corrected for the 
track finding efficiency and have an error (statistical 
plus systematic) of around 5%. 



per unit pseudo-rapidity with px > 0.5 GeV (multiply 
by 3 to get 360°, divide by 2 for the two units of 
pseudo-rapidity, multiply by 1.09 to correct for the 
track finding efficiency). We know that if we include 
all pT that there are roughly 4 charged particles per 
unit rapidity in a "soft" pp collision at 1.8 TeV, and 
the data show that in the underlying event of a "hard 
scattering" there are about 3.8 charged particles per 
unit rapidity in the region pT > 0.5 GeV! If one 
includes all pT values then the underlying event has 
a charge particle density that is at least a factor 
of two larger than the 4 charged particles per unit 
rapidity seen in "soft" pp collisions at this energy. As 
can be seen in Fig. ^ the charged particle density 
in the "transverse" region is a function of pT(jet#l) 
and rises very rapidity at low ^^(jet^l) values. The 
"transverse" charged multiplicity doubles in going from 
PT(jet#l) = 1.5 GeV to pT(jet#l) 2.5 GeV and 
then forms an approximately constant "plateau" for 
PT(jet#l) > 6 GeV. 

Fig. ^ and Fig. compare the "transverse" (Nchg) 
and the "transverse" {pTsum)^ respectively, with the 
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Figure 33. Data from Fig. ^ on the average number 
of charged particles (pr > 0.5 GeV and |?7| < 1) as 
a function of pT(jet#l) (leading charged jet) for the 
"transverse" region defined in Fig. ^ compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of HERWIG 5.9, ISAJET 7.32, and PYTHIA 6.115. 
The errors on the (uncorrected) data include both sta- 
tistical and correlated systematic uncertainties. The 
theory curves are corrected for the track finding effi- 
ciency and have an error (statistical plus systematic) of 
around 5%. 
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Figure 34. Data from Fig. ^ on the average scalar 
Pt sum of charged particles {pt > 0.5 GeV and 
I77I < 1) as a function of pT(jet#l) [leading charged 
jet) for the "transverse" region defined in Fig. 
compared with the QCD "hard scattering" Monte- 
Carlo predictions of HERWIG 5.9, ISAJET 7.32, 
and PYTHIA 6.115. The errors on the [uncorrected) 
data include both statistical and correlated systematic 
uncertainties. The theory curves are corrected for the 
track finding efficiency and have an error [statistical 
plus systematic) of around 5%. 



QCD "hard scattering" Monte-Carlo predictions of 
HERWIG 5.9, ISAJET 7.32, and PYTHIA 6.115. 
Fig. |3| and Fig. |36| compare the "transverse" {Nchg) 
and the "transverse" (PTsum) ^ respectively, with three 
versions of PYTHIA(6.115, 6.125, and no multiple 
scattering, see Table ^0|). PYTHIA with no multiple 
parton scattering does not have enough activity in the 
underlying event. PYTHIA 6.115 fits the "transverse" 
(Nchg) the best, but overshoots slightly the "toward" 
{Nchg) in Fig. |2|. ISAJET has a lot of activity in the 
underlying event, but gives the wrong pT(jet#l) de- 
pendence. Instead of a "plateau," ISAJET predicts a 
rising "transverse" (Nchg) and gives too much activity 
at large pT(jet#l) values. HERWIG does not have 
enough "transverse" (PTsum)- 
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Figure 35. Data from Fig. ^ on the average number 
of charged particles (pr > 0.5 GeV and \ri\ < 1) as 
a function of pT(jet#l) [leading charged jet) for the 
"transverse" region defined in Fig. |2j compared with 
the QCD "hard scattering" Monte-Carlo predictions of 
PYTHIA 6.115, PYTHIA 6.125, and PYTHIA with 
no multiple parton scattering (No MS). The errors 
on the (uncorrected) data include both statistical and 
correlated systematic uncertainties. The theory curves 
are corrected for the track finding efficiency and have 
an error (statistical plus systematic) of around 5%. 



We expect the "transverse" region to be composed 
predominately from particles that arise from the break- 
up of the beam and target and from initial-state radi- 
ation. This is clearly the case as can be seen in Fig. ^ 
where the predictions of ISAJET for the "transverse" 
region are divided into three categories: beam-beam 
remnants, initial- state radiation, and outgoing jets 
plus final-state radiation. It is interesting to see that 
it is the beam-beam remnants that arc producing the 
approximately constant "plateau" . The contributions 
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Figure 36. Data from Fig. |26| on the average scalar pT 
sum of charged particles {px > 0.5 GeV and |ry| < 1) 
as a function of pT(jet#l) (leading charged jet) for the 
"transverse" region defined in Fig. ^ compared with 
the QCD "hard scattering" Monte-Carlo predictions of 
PYTHIA 6.115, PYTHIA 6.125, and PYTHIA with 
no multiple parton scattering (No MS). The errors 
on the (uncorrected) data include both statistical and 
correlated systematic uncertainties. The theory curves 
are corrected for the track finding efficiency and have 
an error {statistical plus systematic) of around 5%. 



from initial-state radiation and from the outgoing hard 
scattering jets both increase aspT(jet#l) increases. In 
fact, for ISA JET it is the sharp rise in the initial-state 
radiation component that is causing the disagreement 
with the data for pT(jet#l) > 20 GeV. 

As we explained in Section B, for HERWIG and 
PYTHIA it makes no sense to distinguish between 
particles that arise from initial-state radiation from 
those that arise from final-state radiation, but one can 
separate the "hard scattering component" from the 
beam-beam remnants. For PYTHIA the beam-beam 
remnants include contributions from multiple parton 
scattering as illustrated in Fig. |l^. Fig. ^ and Fig. ^ 
compare the "transverse" {Nchg) with the QCD "hard 
scattering" Monte-Carlo predictions of HERWIG 5.9 
and PYTHIA 6.115, respectively. Here the predictions 
are divided into two categories: charged particles that 
arise from the break-up of the beam and target {beam- 
beam remnants), and charged particles that result from 
the outgoing jets plus initial and final-state radiation 
{hard scattering component). As was the case with 
ISAJET the beam-beam remnants form the approx- 
imately constant "plateau" and the hard scattering 
component increase as pT(jet#l) increases. However, 
the hard scattering component of HERWIG and 
PYTHIA does not rise nearly as fast as the hard 
scattering component of ISAJET. This can be seen 



Figure 37. Data from Fig. ^ on the average number 
of charged particles {pr > 0.5 GeV and I77I < 1) as 
a function of pT(jet#l) {leading charged jet) for the 
"transverse" region defined in Fig. |2j compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of ISAJET 7.32. The predictions of ISAJETare 
divided into three categories: charged particles that 
arise from the break-up of the beam and target 
{beam-beam remnants), charged particles that arise 
from initial-state radiation, and charged particles that 
result from the outgoing jets plus final-state radiation 
(see Fig. |l^). The errors on the (uncorrected) data 
include both statistical and correlated systematic un- 
certainties. The theory curves are corrected for the 
track finding efficiency and have an error {statistical 
plus systematic) of around 5%. 



clearly in Fig. |40| where we compare directly the hard 
scattering component {outgoing jets plus initial and 
final-state radiation) of the "transverse" (Nchg) from 
ISAJET 7.32, HERWIG 5.9, and PYTHIA 6.115. 
PYTHIA and HERWIG are similar and rise gently 
as pT(jet#l) increases, whereas ISAJETproduces a 
much sharper increase as pT(jet#l) increases. 

There are two reasons why the hard scattering 
component of ISAJET is different from HERWIG and 
PYTHIA. The first is due to different fragmentation 
schemes. ISAJET uses independent fragmentation, 
which produces too many soft hadrons when partons 
begin to overlap. The second difference arises from the 
way the QCD Monte-Carlo produce "parton showers" . 
ISAJET uses a leading-log picture in which the par- 
tons within the shower are ordered according to their 
invariant mass. Kinematics requires that the invariant 
mass of daughter partons be less than the invariant 
mass of the parent. HERWIG and PYTHIA modify 
the leading-log picture to include "color coherence 
effects" which leads to "angle ordering" within the 
parton shower. Angle ordering produces less high px 
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Figure 38. Data from Fig. ^ on the average number 
of charged particles {pr > 0.5 GeV and \r]\ < 1) as 
a function of pT(jet#l) {leading charged jet) for the 
"transverse" region defined in Fig. ^ compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of HERWIG 5.9. The predictions of HERWIG are 
divided into two categories: charged particles that arise 
from the break-up of the beam and target [beam-beam 
remnants), and charged particles that result from the 
outgoing jets plus initial and final-state radiation [hard 
scattering component) (see Fig. ^9|). The errors on 
the (uncorrected) data include both statistical and 
correlated systematic uncertainties. The theory curves 
are corrected for the track finding efficiency and have 
an error [statistical plus systematic) of around 5%. 



radiation within a parton shower which is what is seen 
in Fig. Without further study, we do not know 
how much of the difference seen in Fig. 1^ is due to 
the different fragmentation schemes and how much is 
due to the color coherence effects. 

The beam-beam remnant contribution to the "trans- 
verse" (Nchg) is different for each of the QCD Monte- 
Carlo models. This can be seen in Fig. ^ where we 
compare directly the beam-beam remnant component 
of the "transverse" (Nehg) from ISAJET 7.32, HER- 
WIG 5.9, PYTHIA 6.115, and PYTHIA with no 
multiple parton interactions. Since we are considering 
only charged particles with py > 0.5 GeV, the height 
of the "plateaus" in Fig. ^ is related to the pr dis- 
tribution of the beam-beam remnant contributions. A 
steeper pT distribution means less particles with px > 
0.5 GeV. PYTHIA uses multiple parton scattering to 
enhance the underlying event and we have included 
these contributions in the beam-beam remnants. For 
PYTHIA the height of the "plateau" in Fig. |l| 
can be adjusted by adjusting the amount of multiple 
parton scattering. HERWIG and ISAJET do not 
include multiple parton scattering. For HERWIG and 



Figure 39. Data from Fig. ^ on the average number 
of charged particles [pr > 0.5 GeV and I77I < 1) as 
a function of pT(jet#l) [leading charged jet) for the 
"transverse" region defined in Fig. |2j compared with 
the QCD "hard scattering" Monte-Carlo predictions 
of PYTHIA 6.115. The predictions of PYTHIA are 
divided into two categories: charged particles that 
arise from the break-up of the beam and target 
[beam-beam remnants), and charged particles that 
result from the outgoing jets plus initial and final-state 
radiation [hard scattering component). For PYTHIA, 
the beam-beam remnants include contributions from 
multiple parton scattering (see Fig. |l^). The errors 
on the [uncorrected) data include both statistical and 
correlated systematic uncertainties. The theory curves 
are corrected for the track finding efficiency and have 
an error [statistical plus systematic) of around 5%. 



ISAJET the height of the "plateau" can be adjusted 
by changing the pt distribution of the beam-beam 
remnants. 

8.6. Summary and Conclusions 

For pp collisions at 1.8 TeV min-bias does not 
necessarily imply "soft" physics. There is a lot of 
QCD "hard scattering" in the Min-Bias data. We 
have studied both "local" leading jet observables and 
"global" observables where to fit the data the QCD 
Monte-Carlo models have to correctly describe the 
entire event structure. Our summary and conclusions 
are as follows. 

The Evolution of Charge Particle Jets 
Charged particle jets are "born" somewhere around 
prQet) of about 2 GeV with, on the average, about 2 
charged particles and grow to, on the average, about 
10 charged particles at 50 GeV. The QCD "hard 
scattering" models describe quite well [although not 
perfectly) "local" leading jet observables such as the 
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Figure 40. QCD "hard scattering" Monte-Carlo 
predictions from HERWIG 5.9, ISAJET 7.32, and 
PYTHIA 6.115 of the average number of charged 
particles (pT > 0.5 GeV and \r]\ < 1) as a function 
ofpT(jet#l) {leading charged jet) for the "transverse" 
region defined in Fig. |2j arising from the outgoing jets 
plus initial and finial-state radiation (hard scattering 
component). The curves are corrected for the track 
finding efficiency and have an error {statistical plus 
systematic) of around 5%. 
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Figure 41. QCD "hard scattering" Monte- 

Carlo predictions from HERWIG 5.9, ISAJET 7.32, 
PYTHIA 6.115, and PYTHIA with no multiple 
parton scattering (No MS) of the average number 
of charged particles (pr > 0.5 GeV and \r\\ < 1) 
as a function of p^Qet^l) {leading charged jet) for 
the "transverse" region defined in Fig. 7 arising from 
the break-up of the beam and target {beam-beam 
remnants). For PYTHIA the beam-beam remnants 
include contributions from multiple parton scattering 
(see Fig. |l^). The curves are corrected for the track 
finding efficiency and have an error {statistical plus 
systematic) of around 5%. 



multiplicity distribution of charged particles within 
the leading jet, the "size" of the leading jet, the 
radial flow of charged particles and pt around the 
leading jet direction, and the momentum distribution 
of charged particles within the leading jet. In fact, 
the QCD "hard" scattering Monte-Carlo models agree 
as well with 2 GeV charged particle jets as they do 
with 50 GeV charged particle jets! The charge particle 
jets in the Min-Bias data are simply the extrapolation 
{down to small pr) of the high px jets observed in the 
JET20 data. For a fixed pT(hard), the QCD "hard" 
scattering cross section grows with increasing collider 
energy. As the center-of-mass energy of a pp collision 
grows, "hard" scattering becomes a larger and larger 
fraction of the total inelastic cross section. At 1.8 TeV 
"hard scattering" makes up a sizable part of the "hard 
core" inelastic cross section and a lot of min-bias events 
have 2 TeV or 3 GeV jets. 

The "Underlying Event" 

A hard scattering collider event consists of large 
Pt outgoing hadrons that originate from the large pt 
partons {outgoing jets) and also hadrons that originate 
from the break-up of the proton and antiproton {beam- 
beam remnants). The "underlying event" is formed 
from the beam-beam remnants, initial-state radiation, 
and possibly from multiple parton interactions. Our 
data show that the charged particle multiplicity and 
scalar pT sum in the "underlying event" grows very 
rapidly with the transverse momentum of the leading 
charged particle jet and then forms an approximately 
constant "plateau" forpy(jet#l) > 6 GeV. The height 
of this "plateau" is at least twice that observed in "soft" 
collisions at the same corresponding energy. 

None of the QCD Monte-Carlo models we examined 
correctly describe all the properties of the under- 
lying event seen in the data. HERWIG 5.9 and 
PYTHIA 6.125 do not have enough activity in the 
underlying event. PYTHIA 6.115 has about the right 
amount of activity in the underlying event, but as a 
result produces too much overall charged multiplicity. 
ISAJET 7.32 has a lot of activity in the underlying 
event, but with the wrong dependence on ^^(jet^l). 
Because ISAJET uses independent fragmentation and 
HERWIG and PYTHIA do not, there are clear 
differences in the hard scattering component {mostly 
initial-state radiation) of the underlying event between 
ISAJET and the other two Monte-Carlo models. Here 
the data strongly favor HERWIG and PYTHIA over 
ISAJET. 

The beam-beam remnant component of both 
ISAJET 7.32 and HERWIG 5.9 has the wrong pT 
dependence. ISAJET and HERWIG both predict 
too steep of a pr distribution. PYTHIA does a 
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better job, but is still slightly too steep. It is, of 
course, understandable that the Monte-Carlo models 
might be somewhat off on the parameterization of 
the beam-beam remnants. This component cannot 
be calculated from perturbation theory and must be 
determined from data. With what we have learned 
from the data presented here, the beam-beam remnant 
component of the QCD "hard scattering" Monte-Carlo 
models can be tuned to better describe the overall 
event in pp collisions. 
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